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THE DEVELOPMENT OF A STOCHASTIC MODEL OF THE ATMOSPHERE BETWEEN 

30 AND 90 KM TO BE USED IN DETERMINING THE EFFECT OF ATMOSPHERIC 
VARIABILITY ON SPACE SHUTTLE ENTRY PARAMETERS 

by 

Janet W. Campbell 
(ABSTRACT) 

A stochastic model of the atmosphere between 30 and 90 km was de- 
veloped for use in Monte Carlo Space Shuttle entry studies* The model 
is actually a family of models, one for each latitude- seas on category 
as defined in the 1966 U. S. Standard Atmosphere Supplements. Each 
latitude- seas on model generates a pseudo-random temperature profile 
whose mean is the appropriate temperature profile from the Standard 
Atmosphere Supplements. The standard deviation of temperature at each 
altitude for a given latitude- season model is estimated from sounding- 
rocket data. Departures from the mean temperature at each altitude are 
produced by assuming a linear regression of temperature on the solar 
heating rate of ozone. A profile of random ozone concentrations is 
first generated using an auxiliary stochastic ozone model, also develop- 
ed as a part of this study, and then solar heating rates are computed 
for the random ozone concentrations. 

Pressure and density profiles are calculated from each temperature 
profile by solving simultaneously the hydrostatic equation and the 
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equation of state. Isopynic layers, regions where atmospheric density 
is almost constant, are created at 2h and 90 km altitudes by using these 
two altitudes as boundary points in integrating the hydrostatic equation. 

A computer tape of hb2 sounding-rocket measurements of the atmos- 
phere above 35 km was used to estimate parameters in the temperature 
distribution, and to compare the data's statistical characteristics with 
those of the model. The rocket soundings were sorted according to lat- 
itude-season categories and adjustments were made to assure independence 
of profiles. Because the resultant sample sizes were small, confidence 
intervals associated with the sample statistics were so wide that com- 
parisons of these statistics with model statistics were inconclusive. 

The model was used in Monte Carlo simulations of Shuttle entries to 
study the effect of atmospheric variability on the Shuttle entry para- 
meters : maximum dynamic pressure, maximum stagnation-point heating 

rate, maximum g-load, final down-range distance, and final cross-range 
distance. A sample of 1000 entries was generated for each of the four 
seasons in the 30° latitude band, and the resultant parameter distri- 
butions were analyzed to determine parameter values. In general, autumn 
and winter parameter distributions showed more dispersion than spring 
and summer distributions , because the autumn and winter atmosphere models 
were more variable than the spring and summer models. Pearson distri- 
butions were fitted to the empirical distributions, and design values 
obtained in this manner were compared with the traditional "three- sigma" 
design values. 


t 

I 



ACKNOWLEDGMENTS 


For their approval and support of her graduate program, the author 
wishes to thank the management staff of the National Aeronautics and 
Space Administration’s Langley Research Center, and, in particular, the 
Graduate Studies Committee, the Training Office, and her supervisors in 

the Space Systems Division. 

Individuals who deserve specific acknowledgments are: 

•Dr. Richard G. Krutchkoff, the author's major advisor, for his 
assistance with the stochastic modeling and simulation techniques, 

• Mr. Paul F. Holloway, NASA - Langley, for suggesting the topic and 
for his assistance with the Shuttle trajectory portion, 

•Dr. Walter B. Olstad, NASA - Langley, for his review of the disser- 
tation, particularly the material relating to atmospheric physics, 

•Dr. 0. E. Smith and his associates at NASA’s Marshall Space Flight 
Center for their advice concerning atmospheric modeling, 

• Dr. James M. Russell, NASA - Langley, for suggesting the use of 

Green's ozone formula, 

•Dr. Robert E. Boughner, NASA - Langley, for his assistance with 

the radiative heating model, and 

•Dr. James F. Campbell, the author's husband, for believing in 

Women's Liberation. 


ii 



TABLE OF CONTENTS 


Page 

Title Page . « . ^ 

Acknowledgments • • • H 

Table of Contents 

Chapter 

I. Introduction 1 

II. Atmospheric Models and Modeling Techniques 7 

A general description of the atmosphere ... 7 

Physical laws relating atmospheric temperature, 

pressure and density ...... 12 

Deterministic empirical models . 13 

Deterministic theoretical models » • l6 

Stochastic models 19 

III. The Atmospheric Temperature Model ..... 25 

The source of temperature variation .... 25 

Radiative heating rates * 30 

Model linking temperatures to solar heating rates . 33 

IV. The Ozone Model ...... 39 

The ozone model below 50 km ^0 

The ozone model above 50 km U8 

V. Comparison of the Model with Data 59 

The MRN data tape 59 

Comparison of model temperatures with data 6l 

Comparison of model densities and pressures with data .... 7*+ 

iii 



iv 


TABLE OF CONTENTS (continued) 


Chapter 

Conclusions .... 

VI. The Effect of Atmospheric Variability on Five Important 

Space Shuttle Entry Parameters: An Application of the 

Model 

A two-dimensional random density ‘model 

Procedure . . 

Results and discussion ... 

Conclusions ....... 

VII. Summary and Conclusions 

Summary * 

Concluding observations 


Page 
. 83 


. 85 

. 86 
. 90 

. 92 

. 101 + 
. 106 
. 106 
. 110 


References * 

Appendix A. Equations for Calculating Atmospheric Properties 

Governing equations . 

Solving for p and P given a temperature- altitude 
profile . . 

Solving for z and p given a temperature-pressure 
profile ...... 

Appendix B. Equations for Calculating Solar Heating Rates . 

Ozone absorption coefficients . . 

Solar heating rates 

Calculating optical thicknesses 

Appendix C. Estimating the Parameters in Green's Formula 

Using Ozonagrams 


, . 113 
. . 117 
. . 117 

. . 118 

. . 120 
. . 122 
. . 123 
. . 125 
. . 127 

. . 130 


Criteria for a good fit 


131 



V 


TABLE OF CONTENTS (concluded) 

Page 

Calculating the shaping parameter A 131 

Calculating PO* and P # 132 

Appendix D. Equations for Fitting Pearson Distri hut ions to 

Data * * 135 

The Pearson Type I distribution 139 

The Pearson Type II distribution 1^0 

The Pearson Type III distribution . 1^1 

The Pearson Type IV distribution 1^1 

The Pearson Type V distribution 1^2 

The Pearson Type VI distribution 1^3 

The Pearson Type VII distribution 1^ 

The Normal distribution * ll^ 



I. INTRODUCTION 


In evaluating the performance of various aerospace entry vehicle 
designs, a major tool has been the use of entry trajectory computer 
programs. The majority of these programs use the U. S. Standard Atmos- 
phere, [l], to furnish atmospheric temperatures, densities, and pres- 
sures as a function of altitude. These tables, published in 1962 by 
the Committee on Extension to the Standard Atmosphere (COESA) represent- 
ing 29 U. S. scientific and engineering organizations, were developed 
as a means of standardizing aircraft and spacecraft engineering calcu- 
lations. The basis of the tables is an annual mean temperature profile 
located at 45 ° N. latitude. Thus, the atmospheric properties furnished 
by the U. S. Standard Atmosphere are reasonably close to conditions 
during the spring and fall months at 45 ° N. latitude but are not repre- 
sentative of winter and summer properties nor those at different lati- 
tudes. In 1966, COESA published the U. S. Standard Atmosphere Supple- 
ments, [2], which contains seasonal atmospheric tables for 30°, 45 °, 
60°, and T5° N. latitudes and a set of annual mean profiles for 15° N. 
latitude. 

Although these standard atmospheric tables are valuable for calcu- 
lating nominal trajectories and for providing a basis of comparison for 
engineering calculations made by different individuals, they are inade- 
quate for determining the impact of extreme atmospheres on a vehicle’s 
performance. In designing a spacecraft it is necessary to establish 
a performance envelope within which all flight parameters are expected 
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to lie with high probability regardless of environmental extremes. Once 
this envelope is determined , one can establish a set of design criteria 
stating the extreme loads or stresses which an entry vehicle must be able 
to withstand. 

In past design studies on atmospheric entry vehicles, a common 
method for obtaining a performance envelope (e.g., [3] and [4]), was to 
calculate trajectories using both maximum and minimum atmospheric density 
profiles as shown in Figure 1.1. Generally these are profiles in which 
atmospheric density is three standard deviations above or below its mean 
at all altitudes simultaneously. This method has two major disadvantages. 
One is that since density profiles such as these never occur in any real- 
atmosphere, a design parameter based on this method might be overly con- 
servative and, thereby, require unnecessary expense. A second disadvan- 
tage, more critical than the first, is that these extreme density pro- 
files do not produce extremes in all entry parameters. For example, with 
some spacecraft a more severe total heat load is produced when atmospher- 
ic density is extremely low during the initial phase of entry and sudden- 
ly becomes extremely high at lower altitudes [5]* The reason for this 
severity is that the initial low density causes less deceleration than 
is normal, and thus the spacecraft encounters an extremely dense atmos- 
phere while traveling at an unusually high velocity. 

To account for the fact that extremes in the various entry parameters 
are produced by different atmospheric situations, an alternate determin- 
istic approach, [6], has been to determine analytically for each perfor- 
mance parameter the atmospheric profile which produces an extreme in 
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Percent departure from the 1962 Standard density 


Figure 1.1 - "Extreme 11 atmospheres used in entry vehicle design studies. 
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that parameter, and then design the vehicle to withstand that extreme. 

As with any deterministic approach, however, this has the disadvantage 
that any specific atmospheric profile has a zero probability of occur- 
rence, and thus the design may be overly conservative. Furthermore, 
the degree of conservatism cannot be ascertained since no knowledge is 
provided as to the probability of encountering atmospheric profiles 
similar "to "the design atmosphere. 

The disadvantages associated with deterministic methods are strong 
justification for using statistical methods in establishing design cri- 
teria. A Monte Carlo entry study based on realistic random atmospheres 
can be used to estimate the statistical distribution of any entry para- 
meter. Then a parameter's design value can be selected for any desired 
risk or exceedance probability. The optimum design criterion is one in 
which the risk associated with its design value is a compromise between 
safety and cost. That is, its risk should be low enough to provide a 
high measure of safety but not so low as to make the design impossible 
or unnecessarily expensive. 

The next major space transportation system of the future, the 
Space Shuttle, is currently in the design stage. One major innovation 
in this new system is that the Shuttle will be capable of returning from 
orbit and landing much like today’s conventional aircraft. Thus, the 
Shuttle can be reused, and the waste of disposing with spacecraft will 
be eliminated. With the advent of this reusable spacecraft, however, 
it has become much more important to have a good estimate of the proba- 
bility of exceeding any design value. For example, when a spacecraft 
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was used only once it was not so critical to distinguish between a fail- 
ure probability of .005 and one of .001. With the Shuttle having an ex- 
pected lifetime of 100 missions, however, the difference between single- 
flight failure probabilities of .005 and .001 makes the difference between 
a 40 percent and a 10 percent chance of at least one failure during the 
lifetime of the Shuttle. This, therefore, is further justification for 
using the best available statistical techniques in design studies to ac- 
quire as much confidence as possible in reliability estimates. 

In order to do Monte Carlo simulations of Shuttle entries , it was 
necessary to develop a stochastic atmosphere model capable of rapidly 
generating realistic atmospheric profiles. The author has accomplished 
this for the altitude region between 30 and 90 km, the region of great- 
est impact on Shuttle entry design considerations. This model uses a 
combination of empirical and theoretical techniques to derive a pseudo- 
random altitude profile of atmospheric temperature. The mean tempera- 
ture at each altitude is defined empirically, and variations about the 
mean are explained, in part, by theoretical considerations. An auxil- 
iary stochastic model of atmospheric ozone concentrations, also develop- 
ed as a part of this study, is used to produce variations in solar 
heating rates, which are then linked to temperature variations via a 
linear regression model. Since the three thermodynamic properties of 
the atmosphere, temperature, density, and pressure, are related by two 
physical equations, it is possible to calculate any two quantities given 
an altitude profile of the third. Therefore, this model is fundamental- 
stochastic model of atmospheric temperature, and pressure and 
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density profiles are obtained by using the pseudo-random temperature 
profiles provided by the model to solve tvo deterministic equations. 

A review of atmospheric modeling techniques and current models is 
presented in Chapter II. Chapter III describes the development of the 
temperature model, and Chapter IV presents the auxiliary stochastic 
ozone model mentioned above. Chapter V presents the statistical char- 
acteristics of all three atmospheric profiles generated by the model 
(i.e., temperature, density, and pressure), and compares these with data 
obtained from Meteorological Rocket Network soundings. Chapter VI 
describes a Monte Carlo Shuttle entry study based on the present model. 
The model was used to furnish one pseudo-random density profile at the 
beginning of each entry simulation. In order to be realistic, since 
the Shuttle follows a long shallow entry path, it was necessary to add 
a horizontal variation to atmospheric density. TMb is also discussed 
in Chapter VI. Chapter VII begins with a summary of what has been 
accomplished, points out directions for improvement and future work, 
and states several conclusions regarding the overall approach. 



IX. ATMOSPHERIC MODELS AND MODELING TECHNIQUES 
SURVEY OF LITERATURE 

A general description of the atmosphere 

For the purposes of entry studies it will he assumed that the 
atmosphere extends to an altitude of 90 km (~300,000 ft). Aside from 
the fact that only .0001 percent of the atmosphere lies outside this 
radius, 90 km is a convenient upper limit for the atmosphere because at 
that altitude the relative proportions of the various constituents begin 
to vary, making it more complicated to relate temperature, density, and 
pressure. Below 90 km the molecular weight of air is essentially con- 
stant (28.964 kg), reflecting a constancy in the relative proportions 
of its principal constituents (see Table 2.1). As the molecular weight 
of air begins to vary above 90 km, this quantity becomes a fourth vari- 
able in the equation of state. Furthermore, the hydrostatic equation 
relating pressure and density begins to fail, and, therefore, the math- 
ematical relationships between temperature, density and pressure used 
below 90 km become invalid above that altitude. 

Figure 2.1 shows a general schematic diagram of the atmosphere and 
its various layers. The temperature profile shown here is that which 
forms the basis of the U* S. Standard Atmosphere, [l]. It is a piece- 
wise linear approximation to a mean temperature profile obtained by 
averaging over temperatures from all seasons and latitudes . The various 
layers of the atmosphere indicated in Figure 2.1 are distinguished 
according to whether the rate of change of temperature with height is 
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Table 2.1 - Composition of Dry Air 
Belov 90 km 


Constituent gas and 
formula 

Content , 
percent by 
volume 

— 

Content variable 

relative to its 

normal 

Molecular 

weight* 

Nitrogen (Ng) 

78.084 

- 

28.0134 

Oxygen (Og) 

20. 9 W 

- 

31.9988 

Argon (Ar ) 

0.93b 

- 

39.948 

Carbon dioxide (COg) 

0 . 31 U 

t 

44.00995 

Neon (He) 

0.001818 

- 

20.183 

Helium (He) 

0.000524 

- 

4.0026 

Krypton (Kr) 

0.000114 

- 

83.80 

Xenon (Xe) 

0.0000087 

- 

131.30 

Hydrogen (Hg) 

0.00005 

? 

2.01594 

Methane (CH^) 

0.0002 

t 

16.04303 

Nitrous oxide (NgO) 

0.0005 

- 

44.0128 

Ozone (CO 

0 to .00002 

t 

47.9982 

Sulfur dioxide (SOg) 

0 to 0.0001 

t 

64.0628 

Nitrogen dioxide (NOg) 

0 to 0.000002 

■ 

46.0055 

Ammonia (NH^) 

0 to trace 


17.03061 

Carbon monoxide (CO) 

0 to trace 

+ 

28.01055 

Iodine (Ig) 

0 to 0.000001 

t 

253.8088 


On basis of carbon-12 isotope scale for which C =12. 

, ^The content of the gases marked with a dagger may undergo signifi- 
cant variations from time to time or from place to place relative to the 
normal indicated for those gases. 
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negative or positive. 

The troposphere is the region beginning at the earth’s surface in 
which temperatures decrease as altitude increases. This region contains 
approximately 99.8 percent of the atmosphere’s precipitable water and, 
accordingly, contains the weather systems which affect surface conditions. 
The stratosphere is the region of increasing temperatures between "10 km 
(the tropopause) and "“50 km (the stratopause ) , and the mesosphere, ex- 
tending to "90 km (the mesopause), is a second region of decreasing 
temperatures. The lowest temperatures in the atmosphere occur at the 
mesopause. Up to an altitude of about 65 km, summer temperatures exceed 
winter temperatures. However, between 65 and 110 km summer temperatures 
are lower than winter temperatures with an annual difference of as much 
as 60° K occurring at the mesopause. 

The isother mal tropopause, stratopause and mesopause regions de- 
picted on Figure 2.1 do not represent the shape of real temperature 
profiles. They result from the averaging of temperature profiles in 
which the tropopause, stratopause and mesopause occur at distinct varia- 
ble altitudes. 

The thermosphere begins at the mesopause and extends upward to an 
undefined height. Here temperatures increase exponentially with alti- 
tude until they reach an asymptotic value, known as the exospheric 
temperature , which ranges between 600 and 2000° K depending on the 
amo unt of solar activity. 

This partioning of the atmosphere into a troposphere, stratosphere, 
mesosphere, and thermosphere is based on the temperature structure. 



Another classification dividing the atmosphere into an ionosphere and a 
neutrosphere , above and below 80 km, respectively, is based on the fact 
that above ~80 km the air is highly ionized by high-frequency solar 
radiation. Another dual system divides the atmosphere at 90 km accord- 
ing to whether the molecular weight of air is constant (homosphere) or 
variable (heterosphere). In this study, the temperature-based system 
will be used. 

The shape of the temperature profile can be explained in terms of 
the various heating and cooling mechanisms which exist at different 
altitudes. The temperature ’’bulge” in the stratosphere and mesosphere 
results from the absorption of solar radiation by ozone, a gas which is 
formed by photochemical processes in that region of the atmosphere. Al- 
though ozone is a minor constituent, in terms of its relative concentra- 
tion, it plays a vital role in controlling the nature of our atmosphere. 
It absorbs virtually all of the sun's extremely high frequency radiation, 
wavelengths less than ~3000 A , and thereby shields the earth's surface 
from potentially harmful ultraviolet radiation, X-rays, and gamma rays. 

Mechanisms which create the tropospheric temperature structure are 
more complicated than those in the stratosphere and mesosphere. Approx- 
imately 50 percent of all solar radiation reaches the earth's surface, 
and 90 percent of this incident radiation is absorbed. Heat from the 
surface is then transfered to the air adjacent to the surface, and water 
at the surface evaporates. As this moist heated air rises, it cools, 
becomes saturated, and eventually its moisture condenses, thus releasing 
heat into the atmosphere. These vertically moving air masses, together 
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vith horizontal winds and radiative processes, govern the temperatures 
in the troposphere. 

This combination of conductive, convective, and radiative processes, 
plus the heat exchanges resulting from evaporation and condensation, make 
the troposphere an extremely difficult region to model. Fortunately, 
however, the region of most importance for Shuttle entry considerations 
is the region between 30 and 90 km, the upper stratosphere and mesosphere, 
where maximum heating rates, g— loads, and dynamic pressures occur. The 
present discussion of modeling techniques and current models will be 
limited to those which pertain, at least partially, to this region. 

Physical laws relating atmospheric temperature, pressure, and density 

There are two well-established laws relating the thermodynamic 
properties of temperature, T , pressure, P , and density, p . The 
equation of state of an ideal gas (Perfect Gas Law) is given by 



where K is the universal gas constant and M is the molecular weight 
of air. The hydrostatic equation, relating the rate of change in atmos- 
pheric pressure to density, is expressed as 

dP = -gp dz (2.2) 

where z is altitude and g is the acceleration due to gravity, a 
deterministic quantity which depends on the distance from its point of 
application to the center of the earth. The hydrostatic equation 
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expresses the fact that the difference in pressure between two altitudes 
is equal to the weight of air in a vertical column of unit cross-section- 
al area contained between those altitudes. This relationship assumes 
that the atmosphere is static relative to the earth, an assumption which 
may be made below ~90 km. 

Appendix A gives the simultaneous solution of (2.1) and (2.2) which 
is used to derive the pressure and density profiles corresponding to a 
given temperature profile. It also contains an equation used in this 
study for calculating g aB a function of latitude and altitude. 

Most atmospheric models either assume or derive one atmospheric 
profile (e.g. , temperature, density, or pressure) and then calculate the 
remaining two using the equation of state, (2,l), and the hydrostatic 
equation, (2.2). The models vary according to which profile is the one 
used as a basis and according to whether that profile is empirically 
constructed from experimental data or is derived from a set of theoretical 
laws. The two latter classifications, empirical versus theoretical, will 
be discussed separately. 

Deter™ nistic empirical models 

Probably the best known empirical atmosphere models are the 1962 U. 

S. Standard Atmosphere, [l], and the 1966 U. S. Standard Atmosphere 
Supplements, [2]. These are temperature -based models which assume 
piecewise linear temperature profiles approximating the means of exper- 
imental data. These assumed temperature profiles are then used to tabu- 
late values of pressure, density, speed of sound, coefficient of viscosity, 
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thermal conductivity and other gas properties for altitudes at one- 
kilometer intervals between 0 and 120 km. These tables are used for 
general engineering purposes. For more specialized applications such 
as trajectory analyses of launches from Cape Kennedy, Vandenburg, Wallops 
Island, etc. , reference atmospheres for these and other sites have been 

tabulated, (e.g. , [T]» [8], [9])* 

Empirical modelB are sometimes called "statistical" models when, in 
addition to including means of the various properties, they also include 
standard deviations and other sample statistics. An example of such a 
"statistical" model is reference [10] where Theon, et al., present a 
summary of statistics based on 208 temperature soundings made between 

I960 and 1968. 

Many so-called "models" are actually data summaries such as 
reference [11]. This is an example of the monthly data reports on 
Meteorological Rocket Network soundings at selected sites around the 
world. Besides giving the individual profiles measured during the 
month, they also list, for each site, the means and standard deviations 
of all the accumulated data taken during that month since 1961. The 
various Meteorological Rocket Network sites are shown in Figure 2.2. 

As the quantity of data increases each year, global atmosphere 
models are being developed which account not only for seasonal and 
latitudinal variations, but also for longitudinal or time-of-day vari- 
ations. The model presented by Weidner, Chambers and Lou, [12] is a 
preliminary one resulting from a study initiated in 1968 to develop a 
global model of the atmosphere above 25 km for NASA design criteria 
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Figure 2. 2 - Location of Meteorological Rocket Network launch sites. 
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purposes. This con he used with another global model developed by 
Spiegler and Fowler, [13] » for the first 25 km. Another model develop- 
ed specifically to provide a three-dimensional atmosphere for long hori- 
zontal Shuttle trajectories is that of Bowman [lU]. 

Deterministic theoretical models 

The models described above were classified as "empir ical" because 
they rely largely on experimental data to provide values of temperature, 
density and pressure as a function of height. In some cases, [ll], they 
are simple tabulations of data, whereaB in Other cases, [lU], polynom- 
ials or other such curves have been fitted to the data to give mathemat- 
ical representations of the atmospheric profiles. Such empirical models 
are valuable as engineering tools in that they provide the best estimate 
of atmospheric conditions. Theoretical models, on the other hand, such 
as those of Goody, [15]* Harris and Priester, [16], and Kuhn, [17] > 
attempt to explain observed properties using physical laws and hypotheses. 
Although good representations of the atmosphere are achieved, these 
models never completely reproduce the atmosphere's structure, many 
characteristics of which are not yet fully understood. Such models 
are important to the physicist as they enable him to explain observed 
phenomena, but as engineering tools, theoretical models are, in general, 
computationally too cumbersome, and the properties they predict are not 
as accurate as those of empirical models. Only those theoretical 
techniques and models which influenced this development will be discussed 
here. For an excellent general history of atmospheric temperature 
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modeling techniques, the reader is referred to Kuhn, [IT]. 

Probably the most common approach to describing the atmosphere 
theoretically has been to calculate a temperature profile using the 
assumption that the atmosphere is in radiative equilibrium (e.g., 

Manabe and Mailer, [18], and Manabe and Strickler, [19]). Although 
equilibrium temperature profiles of this type approximate observed 
temperature profiles, (see Figure 2.3), they cannot agree perfectly 
because the atmosphere is not actually in radiative equilibrium. In 
reality, there are regions which act as heat sources and sinks thus 
accounting for known global circulation patterns. The atmosphere is, 
however, close to radiative equilibrium, particularly in the strato- 
sphere [15], and thus radiative processes explain, to a great extent, 
the temperature structure. 

Instead of calculating equilibrium temperatures, an alternate 
approach was that used by Kuhn, [17]. He assumed a mean temperature 
profile and then calculated the corresponding infrared heating and 
cooling rates in the upper stratosphere and mesosphere. His Infrared 
heating and cooling rates added to the solar heating rates of Murgatroyd 
and Goody, [20], give a net heat balance for the atmosphere which agrees 

well with known circulation patterns. 

One major characteristic of any modeling technique based on radia- 
tion theory is its computational difficulty. The radiative transfer 
equation is a multiple integral equation requiring a numerical solution. 
Seemingly minor simplifying assumptions can often have disastrous effects 
on the predicted temperature profile or heating and cooling rates, [17]. 
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Standard Atmosphere 

Manabe and Strickler, [19] 

Manabe and Moller , [18] 


Su, [36] 


different boundary conditions 



Temperature , °K 


Figure 2. 3 - Comparison of radiative equilibrium temperature prof i 1 es 
with 1962 Standard temperature profile. 
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Thus, radiative transfer calculations, in general, are not suitable for 
use in a stochastic model required to generate large numbers of temp- 
erature profiles rapidly. However, energy transfer calculations in- 
volving only solar radiation are much simpler, and these have been in- 
corporated into the present stochastic model. 

Stochastic models 

Historically one of the first statistical treatments of atmospheric 
properties was made by Dines, [21], in 1919 » Using the equation of state 
of an ideal gas, (2.1), and the assumption that in the atmosphere 
departures from the mean temperature, density and pressure are small in 
comparison with these means, Dines derived the approximate relations 

C p = C p r (Pp) - C T r(pT) 

C T = C p r(PT) - C p r(pT) (2.3) 

C p = C p r(Pp) + C T r(PT) 

where C , C m , and are the coefficients of variation of p , T , 

p T r 

and P , respectively, and r( • * ) denotes a correlation coefficient 
between the two indicated quantities. 

In 1954, Buell, [22], used these relationships in conjunction with 
the hydrostatic equation to derive an expression for the vertical gradient 
of a(P) , the standard deviation of pressure. His expression is given 
by 

= _ g a(p) r (Pp) 


(2.4) 



20 


where a(p) is the standard deviation of density. According to this 
relationship, a(P) should be maximum or minimum at altitudes where 
r (Pp) = 0. 

Dine’s work has recently been applied by Smith, et al., [23], to 
a modern problem relating to aerospace vehicle design criteria. Smith 
and his associates solved equations (2.3) for the correlation coeffi- 
cients r ( pT ) , r(Pp) , and r(PT)» using experimental measurements of 
C , C m , and CL. at different altitudes. Their results are shown in 
Figure 2.4. They propose that these correlation coefficients be used 
to determine the appropriate combination of P» p, and T when one of 
these properties is taken to be extreme. For example, r(PT) can be used 
to determine the appropriate value of temperature to accompany an ex- 
treme pressure. This approach is valuable in design problems where 
atmospheric properties are needed at a fixed altitude, e.g>, along 
horizontal flight paths. In such cases, vertical atmospheric profiles 
are not known, and hence the hydrostatic equation is not applicable. 

A number of authors have fitted polynomials to various atmospher- 
ic profiles and, by estimating the statistical distributions of the co- 
efficients in their polynomials, they are able to generate random pro- 
files. One example of this is the work of Essenvanger, [24], in which 
he approximates temperature and wind profiles helov 25 km using 
Fourier series and density profiles below 25 km using Tchebycheff 
polynomials. He uses only one random coefficient in each approximation 
and assumes a Weibull distribution for that coefficient. 

Another example of this approach is the recent work of Engler and 
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Figure 2.4 - Coefficients of correlation between atmospheric properties 


as calculated by Smith et al , [23]. 
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Goldschmidt, [25]. Using a set of 6 ? experimentally measured pressure 
profiles, they fitted polynomials of order k = 2 ,..., 6 to profiles of 
log(P) and found that for most profiles a third-order fit was adequate. 
They then determined that the four coefficients, 3 Q * ^ 

3 ^ , used in their third-order polynomials were highly correlated, and 
so they expressed 3 q , 3 -^ and 3 2 as ii near functions of 3^ * 

distribution of 3^ was estimated to be normal. Using their third-order 
polynomial for log(P) , Engler and Goldschmidt were able to express the 
pressure at altitude z as 

P = P 1 (B 3 , z) P 2 (z) (2.5) 

where they call the factor P £ (z) the steady-state pressure, since it 
is independent of 63 , and they called P^Bj , z) the perturbation 
factor. 

As Engler and Goldschmidt point out, one limitation of their model 
is that the 67 profiles used as their data base were for different 
seasons and locations. They speculate that if the data were for a single 
season and location, the correlations between the coefficients 3 Q * 

3 , , and 3^ would increase. However, this author believes that 

just the opposite would be true. That is, high correlations among the 
3 ,’s reflect a systematic variation in the shape of the log(P) profiles 
which may be due to season and latitude effects. If that is the case, 
then in a single season and latitude the variation in log(P) would be 
smaller and more irregular (i.e., less systematic). Hence, the 3^ 3 
would be less correlated. Nevertheless, the necessity of using as many 
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as four random variables instead of one would not be a serious disadvan- 
tage, and this model appears to be quite attractive for Monte Carlo 
simulations. 

Justus and Woodrum, [26], have estimated the statistical distribu- 
tion of "irregular" variations in atmospheric properties. They assumed 
that an atmospheric parameter such as density or pressure can be expres- 
sed as 

F(z) « F (z) + f(z) (2.6) 

o 

where the term F (z) is made up of diumally repeating components 
o 

and steady-state components persisting over a period of several days. 

That is, two values of F (z) measured 2k hours apart would be identical. 

o 

The component f(z) is c all ed the irregular variation, and measurements 
of this were obtained by differencing profiles made at the same time of 
day on consecutive days. Knowing the distribution of f(z) , a random 
F(z) profile could then be generated for any fixed F q (z) . However, 
the variation in f(z) is not sufficient to account for all the variation 
in F(z) even within the same season and location. One would also need 
a model to provide the temporal and special variation in F q (z) . 

A final stochastic model of interest is a multiple regression model 
for atmospheric density between 30 and 110 km developed by Morgenstem 
and Orner, [27]. A stepwise regression procedure was used in selecting 
the optimum set of independent variables which significantly affect 
density. It has been shown by Jacchia, [28], that a number of variables 
have a strong influence on atmospheric density above 200 km. Using the 
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variables in Jacchia's model as a starting point, Morgenstern and Orner 
determined a set of l4 independent variables which have a significant 
correlation with density in the 30 - 110 km region. These are various 
measurements of solar flux, latitude, annual and semi-annual cycles. 

They divided the atmosphere into three quasi-homogeneous regions and 
used a different subset of these lU parameters in their regression equa- 
tion for each region. Their low altitude density model (30 - 50 km) 
required five independent variables; their middle region model (50 - 65 
km) required six, and their high region model (65 - 110 km) required 
two parameters. They claim that their model accounts for 50 percent 
of the variability in density data between 30 and 75 km. 

This model is not complete enough, as given, for generating 
pseudo-random density profiles. It is of interest, however, because, 
unlike the stochastic models discussed above, it is an effort to 
account for atmospheric variation by identifying the major causes of 
that variation. , This was the approach taken in the development of the 
present stochastic temperature model. 



III. THE ATMOSPHERIC TEMPERATURE MODEL 


The single most dominant feature of the altitude region between 30 
and 90 km is the temperature “bulge" caused by the absorption of solar 
energy by ozone. Thus, the temperature profile appears to be a logical 
starting point in attempting to model thiB region. It is for this reason, 
primarily, that the author has chosen to use temperature as the basis of 
the present atmosphere model. Another reason is that it is difficult to 
perturb a standard pressure or density profile and still have the associ- 
ated temperature profile retain its characteristic shape as shown in 
Figure 2.1. 

The source of temperature variation 

To begin modeling variations in atmospheric temperature in the 30 - 
90 km region, the first step is to identify the underlying processes 
which cause that variation and determine which ones are random. Any 
variable which can readily be ascertained without error at the time of 
the spacecraft's entry will be considered nonrandom. Variables which can 
only be estimated statistically at the time of entry will be considered 
random. 

Among the nonrandom causes of temperature variation, a major one 
is nearness to the sun as determined by season and latitude. The present 
model accounts for this source of variation by requiring that season 
and latitude be fixed before estimating the various model parameters. 

That is, the model furnishes pseudo-random temperature profiles 
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characteristic of a specified season and latitude range. 

Another nonrandom source of variation is time of day. r The tempera- 
ture model is based on rates of absorption of solar radiation; therefore, 
it is only representative of profiles during daylight hours. At the 
present time there is insufficient data to estimate model parameters for 
specific times of day. However, the standard deviation of atmospheric 
temperature at each altitude level is specified by the user of the model, 
and thus the model includes all observed variability in daylight tempera- 
tures for a given season and latitude band. 

There are two major sources of temperature variation which are 
considered random. One is the creation and destruction of constituents 

c 

which absorb and/or emit radiation, and the other is the random movement 
o£ air masses. Since the atmosphere is not in radiative equilibrium 
there are regions which act as heat sources (have a positive radiative 
heating rate) and other regions which act as heat sinks (have a negative 
radiative heating rate). Figure 3.1 shows the heating and cooling rates, 
in degrees per day, in the 30 - 100 km region as calculated by Kuhn, [17] • 
These gradients set up global circulation patterns which maintain an 
overall balance in the atmosphere^ heat budget. Superimposed on these 
steady- state circulation patterns are winds or small-scale random air 
movements. A great amount of attention has been focused on random wind 
modeling, (e.g. , [29], [30], [ 31 ] » and [32]). The present temperature 
model does not include the effects of winds explicitly, although there 
is a random error term which might be considered as resulting from winds. 
It may be possible to include winds more directly by combining the present 
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model with the random wind model described in [29]* 

The radiative heating rate at any point is a function of the radiant 
energy being absorbed and emitted at that point. These absorptions and 
emissions depend largely on the amounts of various radiatively active 
gases present. Electromagnetic radiation in the atmosphere has two major 
sources as illustrated in Figure 3«2. Radiation from the sun, with almost 
all wavelengths less than 4 p, lies primarily in the optical spectral 
region, whereas radiation from the earth and its atmosphere, with most 
wavelengths greater than h p, lies in the infrared spectral region. The 
two types of radiation are generally treated separately. 

Solar energy is absorbed primarily by oxygen in the thermosphere, 
ozone in the stratosphere and mesosphere, and water vapor in the trop- 
osphere. Terrestrial radiation (including that from the atmosphere) is 
absorbed by water vapor, carbon dioxide, and ozone. Of these four 
radiatively active gases only water vapor and ozone have concentrations 
which vary significantly from day to day, and of these ozone is the 
major influence in the upper stratosphere and mesosphere. Thus, it was 
decided to model the variations in ozone and use these to produce varia- 
tions in temperature between 30 and 90 km. 

One advantage in linking atmospheric temperatures to radiative 
processes is that radiation, particularly solar radiation, is the only 
obvious means by which temperatures in two widely separated altitude 
regions can be correlated. For example, suppose the ozone in a high 
altitude region absorbs more solar energy than is normal. This means 
that less radiation reaches the ozone below that altitude and, as a 




Figure 3.2 - Spectra] distribution of solar and terrestrial (infrared) 


rad i at i on. 
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consequence, less is absorbed. Since temperatures vary according to 
the amount of energy absorbed, then an increase in temperature (high 
absorption) at a high altitude would tend to be associated with a de- 
crease in temperatures (low absorption) at some lower altitude, and 
vice versa. 

The difficulty in modeling temperature profiles, or any other at- 
mospheric profile, lies in modeling these interlayer correlations. 

Since pressures and densities are obtained by integrating over the 
temperature profile, it is important to generate temperature profiles 
with realistic shapes, i.e., temperature profiles whose deviations from 
the mean in one region are realistic relative to those in other regions. 
Otherwise, it is possible, for example, to construct a temperature model 
having a correct temperature distribution at each individual altitude 
which, at the same time, produces a very unrealistic distribution for 
the integrated temperature profile (i.e., temperature integrated 6ver 
altitude). Since pressure and density depend on this integral, their 
distributions would be adversely affected. 


Radiative heating rates 


A basic premise "underlying the development of the present model 
is that there is a strong positive linear correlation between the temper- 
ature at altitude z , T(z) , and the radiative heating rate at that 
altitude, T(z) . The heating rate, in degrees per unit time , is given 
by 


T(z) 


1 OF 
pc p dz 


( 3 . 1 ) 
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where p is atmospheric density, Cp is the specific heat at constant 
pressure (a constant), and F is the energy flux (the amount of energy 
per unit time passing through a unit area perpendicular to the z di- 
rection) . 

The change in F across a layer of incremental depth dz is 

dF - rate of emission - rate of absorption (3.2) 

The total flux F can be separated into two disjoint spectral 
terms, the flux of optical or solar radiation S and that of infrared 
radiation R . Thus, one can write F=S+R or dF=dS+dR. 

Since the optical energy is basically solar radiation and the atmosphere, 
for most practical purposes, does not emit energy at those wavelengths, 
one can write 

dS = - rate of absorption of solar energy (3.3) 

The solar energy absorbed by a layer depends on the quantity of any 
absorbing gases in that layer (e.g., oxygen, ozone, water vapor) and 
on the amount of direct and reflected solar radiation which is incident 
on that layer. This latter dependency means that dS at altitude z 
depends not only on the quantity of an absorbing gas at altitude z , 
but also on the amounts of that gas at other altitudes, particularly 
at altitudes above z . Thus, in the 30 - 90 km region, variations in 
dS can be determined largely, by modeling variations in the ozone 
profile. 

The change in infrared (IR) flux is given by 
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dfi = rate of IR emission - rate of XR absorption (3**0 

The infrared emission rate depends on both the temperature and pressure 
of the emitting gas as veil as on its Relative concentration. The ab- 
sorption rate depends not only on these properties for the absorbing gas, 
but also on these properties for the surrounding gases. Thus, cLR at 
one altitude is a function of the temperature profile, the pressure 
profile, and the appropriate constituent density profiles. 

For a given ozone profile, a calculation of dS is reasonably 
straight- forward. One reason for this is that absorption coefficients 
needed in the calculation are well-behaved functions of wavelength in 
the visible and ultraviolet spectral regions and are virtually indepen- 
dent of atmospheric temperatures and pressures. An average absorption 
coefficient can be assumed to apply over a fairly wide wavelength band 
without sacrificing much accuracy in the resulting calculations. In 
the infrared region of the spectrum, on the other hand, calculations are 
very difficult. Thousands of overlapping absorption lines, each with 
its own absorption coefficient, are contained in any small spectral 
band. As it is impossible to treat each line separately, a number of 
band models representing various levels of approximation are used, 

(e.g., [33], [3^], and [35]). Besides being mathematically formidable, 
these models involve parameters which are dependent on pressure and 
temperature. As a consequence, calculations involving the infrared 
radiative flux in the atmosphere are highly complex, time consuming, 
and poorly suited for the present purposes. 
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Fortunately, it is possible to concentrate on correlating tempera- 
tures to the solar heating rate 


T s (z) = - 


i as 

pc p dz 


(3.5) 


and ignore the effect of the infrared processes on the variation of 
temperature. Theoretically, such an omission is justifiable provided 
that either dR is insignificant relative to dS, or dB is closely 
correlated to dS so that the partial correlation of T on dR, given 
dS, is small. The former justification is not valid, in general, since 
the atmosphere is close to radiative equilibrium so that dR - -dS. 
However, according to Kuhn's calculations, an increase in ozone produces 
a decrease in dR above UO km and an increase in dR between 30 and 
HO km. Calculations made in the present investigation indicate that the 
same is true of dS. That is, when ozone is increased, dS decreases 
above HO km and increases between 30 and HO km. Thus, it appears that 
changes in dS and dR resulting from ozone variations are 
positively correlated in the regions between 30 and 90 km. 

Calculations of T in the present study are based on a three- 

8 

step non-gray radiative model of the atmosphere used by F. Y. Su, [ 36 ]. 
Details of these calculations are presented in Appendix B. 


Model linking temperatures to solar heating rates 

The decision to relate variations in atmospheric temperature to 
variations in the solar heating rate of ozone required first the 
development of a stochastic model for ozone. Details of this 
development are described in the next chapter. By using this model to 



generate pseudo-random ozone profiles, it is possible to vary the 

heating rates T (z), and thus effect changes in the temperature 
s 

profile T(z). 

To relate T(z) and T (z), a linear regression model is assumed. 

s 

That is, it is assumed that the temperature at altitude z is given 
by a linear function of the general form 

T(z) * 3 Q (z) + B 1 (z)T s (z) + e(z) ( 3 . 6 ) 

where 3 Q (z) and 3 ^(z) are regression coefficients and e(z) is a 
random component representing unmodeled variation (error). 

Equation (3.6) can be written in a more definitive form as 

[T (z) - u (z)3 

T(z) - P T (z) + Xcr T (z) + e(z) (3*7) 

s'* * 

where y T (z) and o T (z) are the mean and standard deviation of T(z), 

respectively; y (z) and o (z) are the mean and standard deviation of 
s s 

, • 

T (z); X is the coefficient of correlation between T(z) and T (z), 
s s 

and e(z) is a normal random number with zero mean and variance given by 

o®(z) = (1 - X 2 )cr 2 (z) ( 3 . 8 ) 

There are no data available which provide simultaneous measurements 

of T(z) and T (z), or of T(z) and ozone, in the 30 - 90 km region. 

5 

Thus, it is not possible to use the usual least-squareB methods for 
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estimating the regression coefficients. Instead, P,p(z) and c^,(z) 

were estimated from sounding- rocket. temperature data. A large volume 

of data is available, [37]* on ozone concentrations in the troposphere 

and stratosphere. These data were utilized to develop a random ozone 

model which was then programmed to generate a large sample of T a (z) 

profiles from which y (z) and o (z) could be estimated. 

s s 

The choice of X = .9 in the present model was based on two consider- 
ations. First, a qualitative comparison was made between real measured 
profiles end model-generated profiles using a range of X values between 
0 and 1, As seen in Figure 3.3* the choice of X determines the amount 
of small scale random scatter in modeled temperature profiles. A value 
of X = 1 eliminates e(z) in equation (3.7) and, therefore, results 
in smooth temperature profiles influenced only by changes in T g . A 
value of X = 0, on the other hand, eliminates the dependence of T 

on T and accounts for all temperature variation with e(z) , thus 

s 

producing a "white noise" type of variation. That is, with X = 0, 
temperatures at any two altitudes would be uncorrelated regardless of 
how close the altitudes are. This comparison of model with data indicat- 
ed that a choice of X between .7 and .9 gave the model an amount of 
s mall -scale random fluctuation resembling that in the data. 

A second consideration in the choice of X was based on the need 
for a sufficiently high correlation between temperatures at adjacent 
altitudes. In the present model pseudo-random temperatures are generated 
at one-kilometer intervals; thus, "adjacent" altitudes are one kilometer 
apart. It follows from equation (3.7) that the correlation between 
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temperatures at any two altitudes cannot exceed X . In the data, 
temperatures at adjacent altitudes generally had very high correlations. 
Thus, to maximize these correlations in the model within the constraint 
,7 < X < .9 , the value X = .9 was selected. Figure 3.^ shows a 
comparison of two measured profiles with two pseudo-random temperature 
profiles using X = .9 . Statistical characteristics of the model are 
compared with those of experimental data in Chapter V. 

Before proceeding to a description of the ozone model, one further 
word of explanation is in order. It will he noticed in the next chapter 
that the ozone model gives ozone concentrations as a function of atmos- 
pheric pressure, P . Quite often atmospheric pressure instead of alti- 
tude is used as a vertical scale. This is possible because P is a 
strictly monotonic (decreasing) function of altitude. In the present 
model, it was found to be more convenient to use P or, more exactly, 
log(P) as the independent variable instead of altitude z . This is 
because the ozone model is based on log(P) . Thus, all the quantities 
dependent on z in equation (3.7) were transformed so as to be functions 
of log(P) so that the model gives a temperature-log(P) profile. Using 
the equation of state and the hydrostatic equation, it is possible to 
solve for z and p as functions of the T-log(P) profile. Details 
are explained in Appendix A. 



Altitude . km Altitude , km 
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IV. THE OZONE MODEL 


Atmospheric ozone has been the subject of much scientific scrutiny 
since its existence was first hypothesized by Schobein in l8Uo, (see 
historical discussion in [38]). Ozone is formed by photochemical proces- 
ses at altitudes where molecular oxygen 0^ end atomic oxygen 0 exist 
simultaneously. Atomic oxygen is formed in the upper atmosphere by the 
dissociation of 0^ when the latter is exposed to high-frequency solar 
radiation. Collisions between 0^ and 0 produce ozone, 0^ . Ozone 
is subsequently destroyed when either 0^ and 0 collide, forming two 
0g molecules, or when 0^ is dissociated into 0 and 0g by solar 
radiation. 

In the mesosphere and upper stratosphere, the rates by which 
ozone is produced and destroyed through these processes balance one 
another so that it should be possible to predict concentrations using 
photochemical equilibrium theory. Unfortunately , however, vyry few 
measurements of ozone have been made in this region, and consequently 
the necessary reaction rates have not been clearly established. 

Beginning at the stratopause and extending downward, photochemical 
equilibrium gradually disappears, and hence theoretical prediction 
methods become less applicable for modeling these ozone concentrations. 
In this region the best predictions of ozone appear to be empirical 
ones based on measured profiles. 

In the present statistical model of ozone, it has been assumed 
that ozone profiles can be described by empirical expressions requiring 
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only a small set of parameters to define any specific profile. It is 
further assumed that the statistical distributions of these parameters 
can he estimated from data. Pseudo-random ozone profiles can then he 
generated hy randomly sampling parameters from their respective distri- 
butions . 

The ozone model below 50 km 

In 1963 a network of observation stations was established by the 
Air Force Cambridge Research Laboratories (AFCRL) to systematically 
measure the vertical distribution of ozone in the troposphere and strat- 
osphere. Figure 4.1 shows the stations making up this "ozonesonde” 
network . 

Figure 4.2 is a typical "ozonagram" published by the AFCRL in a 
four-Volume series entitled "Ozonesonde Measurements over North America”, 
[ 3T ] * The partial pressure of ozone, PO^ , is plotted on the left as 
a function of total atmospheric pressure P . P0^ reaches a maximum 

between approximately 20 and 25 km. Above that point the P0^ profile 
decreases smoothly as it approaches a photochemical equilibrium profile. 
Below its maximum the P0^ profile is characterized by large irregular 
fluctuations in the lower stratosphere and by a constancy or sometimes 
a slight increase in the troposphere. The curve on the right in Figure 
4.2 is a simultaneous reading of atmospheric temperature. Unfortunately 
these curves usually end somewhere between 30 and 35 km, and, thus cannot 
be used for estimating the regression coefficients in equation (3.6) 
where the region of interest is 30 < z < 90 km . 



Station 


Lat.(Deg.N) Long. (Deg. W) 


A I brook Field, Canal Zone 

Colorado State University, Fort Collins 

Eielson AFB, Fairbanks, Alaska 

Florida State University, Tallahassee 

Fort Churchill, Manitoba, Canada 

Goose Bay, Labrador 

L, G. Hanscom Field, Bedford, Mass. 

Thule AFB, Greenland 

University of New Mexico, Albuquerque 

University of Washington, Seattle 

University of Wisconsin, Madison 


9.0 

79.6 

40.6 

105. i 

64.8 

I *>7.9 

30.4 

84.3 

58.8 

94. 1 

53.3 

60.4 

42. 5 

71.3 

76.5 

68.8 

35.0 

106.6 

47.4 

122. 3 

43. 1 

89-4 
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Although ozone concentrations below 25 Km do have some effect on 
solar* heating rates above that altitude, because of reflected solar 
radiation passing upwards through the atmosphere, emphasis in the present 
investigation is placed on modeling the ozone profile above 25 Km. For 
this reason it was not considered necessary to include the large irregu- 


lar variations below the maximum PO^ level. 

An excellent empirical formula was used by Green, [ 39 1 » for 
approximating P0 3 profiles such as the one in Figure 4.2. The formula, 
expressing P0 3 8-s a faction of P , is given by 


PO = 4 PO* -^4- 2 

3 3 (l+e x r 


(U.l) 


where 

x = A log (P/P*) ; (^*2) 


PO^ is the maximum partial pressure of ozone; P* is the total 
atmospheric pressure at the altitude where P0 3 is maximum, and A is 
a shaping parameter. Equation 4.1 will hereafter be referred to as 

Green's formula. Examples of how well Green's formula fits ozonagram 

* 

data are shown in Figure 4,3. In each case the parameters PC> 3 and 
P* were selected so as to match the data above the ozone maximum as 
closely as possible. The four ozonagram shapes are typical of their 
respective stations. In general, both P0 3 and P* increase with 
latitude. 

By integrating over the ozone profile, It is possible to derive 

* 

an approximate expression for total ozone, TQ 3 , as a function of P0 3 



Pressure , mb Pressure . mb 



Altitude , km Altitude , km 










45 


and the shaping parameter A . Setting this expression equal to the 

measured total ozone indicated on each ozonagram, one can solve for A 

# * 

in terms of PO^ and T0 3 . Details as to how PO^ , P* » and A are 

calculated for each ozonagram are contained in Appendix C. 

The statistical distributions of PO^ and P* were estimated from 
the ozonesonde network data in volumes 2 and 3 of reference [37] ■ Each 
ozone network station was assigned to one of five latitude hands, and 
within each latitude hand (except the 15° hand), the ozonagrams were 
divided into four seasonal groups. These season-latitude definitions 
conform to those established by the Meteorological Rocket Network and 
are shown in Table 4.1. This classification results in 17 latitude- 
season categories, and the parameter distributions for each must be es- 
timated. This is a tedious Job, and, to date, has only been accomplished 
for the ozonagrams in the 30° latitude band, the band in which the Shuttle 
landing site (Cape Kennedy) is located. If the need arises it can be 
completed for all 17 categories, and the data in volumes 1 and 4, [37], 
can also be included. 

The parameters PO^ and P* have been estimated for each ozonagram 
using the procedure outlined in Appendix C. Samples of these parameters 
and the TO^ values read directly from the ozonagrams have been tabulat- 
ed for the four season categories in the 30° latitude band, and their 
means, standard deviations and correlations have been estimated. These 
statistics are listed in Table 4.2. 

Since most of the variation in total ozone results from large ir- 
regular fluctuations in the lower stratosphere, variations which are 
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Table 4.1 - Season - Latitude Definitions 


Meteorological Rocket Network Latitude Band Definitions 

15° Band 

0.00° - 22.50° 


30° Band 

22.51° - 37.50° 


45° Band 

37-51° - 52.50° 


60° Band 

52.51° - 67.50° 


75° Band 

67.51° - 90.00° 


Meteorological Rocket Network 

Season Definitions : Northern 

Latitudes 

Spring 

April - May 


Summer 

June - August 


Autumn 

September - October 


Winter 

November - March 


Meteorological Rocket Network 

Season Definitions : Southern 

Latitudes 

Spring 

October - November 


Summer 

December - February 


Autumn 

March - April 


Winter 

May - September 





Ozone Statistics in 30 Latitude Band 


i log<P*} 


Parameter 

Mean 

Standard 

Deviation 


Correlations 




* 

TO 3 

P* TOg tOg(POg) l0g(P*) 


FLORIDA - AUTUMN (sample size = 6) 


PO* 


0,0136 

0.001S19 

1.00 0.13 -0.51 0.14 

3256 

164.3 

0.13 1.00 -0.45 0.19 

0.283 

0.0194 

-0.51 -0.45 1.00 

-4.306 

0.1121 

0.19 1.00 0.20 

8.088 

! 0.0510 

0.14 0.20 1.00 


HEW MEXICO - AUTUMN (sample size * 7) 


0.0132 

0.001463 

1.00 

0.13 

0.39 


0.14 

3661 

164.4 

0.13 

1.00 

0.11 

0.16 


0.259 

0.0119 

0.39 

0.11 

1.00 



-4.329 

0.1102 


0.16 


1.00 

0.17 

7.885 

0.0597 

0.14 



0.17 

1.00 

WINTER (sample size 

* 14) 





0.0155 

0.002146 

1.00 

0.46 

0.60 


0.48 

3151 

419.6 

0.46 

1.00 

0.39 

0.48 


0.303 

0,0286 

0.60 

0.39 

1.00 



-4.175 

0,1426 


0.48 


1.00 

0.49 

8.047 



0.1354 

0.48 



0.49 

1.00 


NEW MEXICO - WINTER (sample size * 22) 


0.153 

0.002519 

1.00 0.48 0.55 0.51 

3349 

812.6 

0.48 1.00 0.8! 0.48 

0.316 

0.0443 

0.55 0.81 1.00 

-4.191 

0.1778 

0.48 1.00 0,52 

8.091 

0.2255 

°- 51 0.52 1.00 



































not accounted for by Green's formula. It was decided to use a fixed 
value for TO^ in our model. Accordingly, TO^ was set equal to its 
latitude-season average and was not allowed to vary. 

The empirical cumulative frequencies of both PO^ end P* were 
compared with normal and log-normal cumulative distribution functions. 

This was done by plotting both the parameter and its natural logarithm 
on normal probability paper, obtaining a least-squares fit, and then 
comparing the sums of squared errors (errors being defined as differ- 
ences between empirical cumulative frequencies and the normal cum u lative 
distribution function). Neither the normal nor the log-normal distribu- 
tion emerged as a clear choice, although the log-normal fit appeared to 

« 

be slightly better in most cases for both PO^ aud P* . Furthermore, 

* 

the log-normal distributions for both PO^ and P* produced smaller, 

more realistic, variances for PO^ at each altitude than did normal 

distributions for either or both parameters. Thus, log-normal distribu- 
te 

tions for PO^ and P* were used. In most dases there was insufficient 
data to warrant fitting a Pearson distribution or attempting a more so- 
phisticated estimate of the parameter distribution. 

The ozone model above 30 km 

Six hundred pseudo-random ozone profiles were generated by sampling 

# 

from the P# and PO^ distributions as described above. Solar heating 

rate profiles T (z) corresponding to each ozone profile were calculated 
s 

* 

and the mean and standard deviation of T were estimated at each alti- 

6 

tude. It was found that the profiles of heating rate averages agreed 





quite well with deterministic T profiles predicted by Su, [ 36 ], 

8 

and Manabe and Strickler, [19 ]. Comparisons are shown in Figure U.h. 

On the other hand. Green's formula and the above parameter distributions 

* 

produced unrealistically large variations in T at altitudes above 

s 

approximately 50 km. At these altitudes, the heating rate is almost 

directly proportional to P0^ * and although the latter is quite small 

* 

in magnitude, the variations in PO^ and P* caused PO^ to vary by 
several orders of magnitude. For example, using Green's formula at an 

1 

altitude of 75 km, the parameter variations result in a range of 
PO^/PO* ratios between 10~ 2 and 10~^. This is sufficient to produce 
a three-order-of-magnitude variation in heating rates at that altitude 
resulting in heating rates as high as 30° per day , a physically un- 
realistic value. 

In reality, ozone in the mesosphere (above 50 km) is in photo- 
chemical equilibrium, for most practical purposes, end its concentrations 
are not as variable as those in the stratosphere. In fact, it is 

it 

unrealistic to allow variations in PO^ and P* to cause appreciable 
changes in PQ^ above 50 km. For this reason, it was decided to mod- 
ify the ozone model in the mesosphere so as to reduce the variation in 
PO^ in that region. 

In Green's formula, (U.l), the term e X approaches 0 as altitude 
increases such that e X « 1 above 50 km. Hence, 

PO = UPO* e x 


( 4 . 3 ) 



Figure 4. 
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T 

s 


prof i les wi th 
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Substituting from equations (U.l) and (^.2), this becomes 

PO^ z exp [A log (P/P*) + B] 

where A is the shaping parameter given in Appendix C as 

A = A PO*/TO_ 

0 0 3 

( for constants A and TO. ) , and 
o 3 

B = log (UPO*) 


(h.k) 


(U.5) 


ik.6) 


The ozone model above 50 km was modified a6 follows; The partial 
pressure of ozone above 50 km is given by 


PO^ = exp [A log (P/P*) + B] 


(U.T) 


where now P* is assumed to be constant and serves strictly as a non- 
dimensionalizing factor. Furthermore, it is now assumed that 


A ~ N (u A , 0 A ) 


(4.8) 


and 


B - N (Ujj * Op 


(U.9) 


Thus, PO^ has a log-normal distribution with mean 


[ 


E(P0 ) = exp log (P/P*) + log 2 (P/P*) + + (^*10) 


2 n 
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and variance 


Var(P0 3 ) = E(P0 3 )' 


j„ P r, 


a 2 log 2 (p/p*) + a 2 


] 


(>*.11) 


The quantity log(P*), which is a normal random number below 50 km, 
is set equal to its expected value above 50 km. 

To maintain continuity with Green's forumula below 50 km, the 
means of A and B are, respectively. 


h A = A q E (F0 3 )/T0 3 


(4.12) 


and 


y. B * log(4) + E [log (PO^)] 


(4.13) 


Since P/P* approaches 0 as altitude increases, log(P/P*) + “» , 


and, therefore, the coefficient of variation of PO^ , 


cv(po 3 ) 


= { eXP \° 


Z k log 2 (P/P*) + cr 2 


H 


1/2 


x 100$ 


(4.l4) 


increases monotonically with altitude. The rate by which this coefficient 

2 2 

of variation increases is controlled by the selection of and . 

2 2 

Accordingly, and were selected so as to satisfy two boundary 

conditions placed on CV (PO^) . At the 50 km boundary the coefficient 
of variation was set equal to the coefficient of variation produced by 
Green's formula, a value which varied between 48$ and 83$ depending on 
the season-latitude model. At the upper boundary, z = 90 km, PO^ was 
assumed to have a coefficient of variation of 100$. These two boundary 
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conditions resulted in a pair of simultaneous equations which were 
2 2 

solved to give a g • 

The four mean PO^ profiles produced "by this model for the 30 
latitude hand are shown in Figure ^.5» Th e dashed lineB are loci of 
maxima and minima occurring in a pseudo-random sample of size 600. Figure 
U.6 compares histograms of PC> 3 values produced hy the model at the 
pressure height of 1000 N/m 2 with corresponding histograms from the data 
in reference [37]. Four mean heating-rate profiles corresponding to the 
P0^ means in Figure U.5 are shown in Figure U.7» again with loci of 
maxima and minima. Table k.3 lists the mean and coefficient of variation 
of heating rate as a function of altitude for these four seasons. 

It should be remembered that this ozone model was developed as a 
tool for modeling temperature variations in the stratosphere and meso- 
sphere. For a large enough X in equation (3.7) » the shape of the 
temperature distribution at any altitude will essentially be that of the 
heating -rate distribution at that altitude. However, the actual magni- 
tudes of the mean and standard deviation of heating rate will not affect 
the mean or standard deviation of temperature, since the two latter 
quantities are built into the model. Thus, further fine-tuning of the o- 
zone model to reduce its variance was not considered necessary. The 
only type of modification which could significantly improve the model 
for its present purposes would be a change affecting the shape of the 
heating-rate distribution; thiB could only be accomplished by changing 
the distribution of PC> 3 . For example, if the distribution of F* is 
estimated to be something other than log-normal, then the PO^ 
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Figure 4.6 - Comparison of model and data histograms of PO 

2 

a pressure height of 1000 N/m . 




Figure 4. 7 " Mean and extreme T profiles based on sample of 600 


model -generated T profiles. 




Table 4,3 - Mean# end Coefficients of Variation T at discrete altitudes 



Spring 

Simmer 

Autumn 

Winter 

r 


Spring 

Sunnier 

f* 

Autumn 

Winter 

Altitude, 

bn 

Mean 

o/day 

CV»S 

Mean 

o/day 

CV.J 

Mean 

o/day 

CV.J 

Mean 

o/day 

, 

cv.j 


Altitude, 

kzn 

Mean 
0 / day 

CV.J 

Mean 

o/day 

CV.J 

Mean 

o/day 

CV.J 

Mean 

o/day 

CV.J 

0 

.ODD 

48.036 

.000 

39.850 

.000 

48.359 

.000 

58.489 

r 

U 6 

4.770 

47.259 

5-338 

37.431 

4.219 

38.740 

4.806 

64.968 











4? 

4. 6U9 

50.571 

5.256 

Uo. 394 

4.069 

41.022 

4.744 

70.067 

1 

.000 

U6. 5^6 

.000 

38.515 

.000 

46.390 

.000 

55-959 


48 

4.511 

53.754 

5.14? 

43.229 

3.910 

44.793 

4.666 

75.139 

2 

.000 

45.045 

.000 

37.156 

.000 

44.384 

.000 

53.421 


49 

4.361 

56.830 

5.019 

45.957 

3.746 

47. 677 

4.574 

. 80.208 

3 

.000 

43.533 

.000 

35.780 

.000 

42.336 

.000 

50 . 880 


50 

3.940 

52.963 

4.603 

42.806 

3-447 

45.886 

3.756 

71.653 

4 

.001 

42.008 

.001 

34.388 

.000 

40.236 

.001 

40.202 











5 

.001 

40.469 

.001 

32.982 

.001 

38.079 

.001 

45.631 


51 

3.741 

54.393 

4.493 

44.326 

3.250 

47.160 

3.552 

73.330 











52 

3.543 

55.696 

4.299 

45.742 

3-059 

46.355 

3.348 

74.799 

6 

.001 

38.917 

.001 

31.562 

.001 

35-862 

.001 

42 , 970 


53 

3.348 

56.896 

4.103 

47.073 

2.873 

49.409 

3.148 

76.069 

7 

.002 

37.350 

.001 

30.127 

-001 

33-577 

.002 

40.206 


54 

3-157 

58.009 

3.907 

48.335 

2.693 

50.574 

2.953 

77.177 

5 

.002 

35.768 

.002 

28.675 

.002 

31.213 

.002 

37.592 


55 

2.972 

59-051 

3.713 

49.538 

2.521 

51.622 

2.764 

78.145 

9 

* DO 3 

34,170 

.003 

27.210 

.002 

28.765 

.003 

34.906 











10 

.005 

32.559 

.004 

25-735 

.003 

26. 227 

.005 

32.254 


56 

2.792 

60.034 

3.522 

50.692 

2.356 

52.639 

2.562 

78.995 











57 

2.620 

60.969 

3-336 

51.807 

2.200 

53.634 

2.408 

79-745 

11 

.007 

30.935 

.006 

24.258 

.005 

23.595 

.007 

29.676 


58 

2.455 

61.866 

3.156 

52.890 

2.051 

54.613 

2.243 

80.413 

12 

a 010 

29.297 

.006 

22.785 

.007 

20.666 

.010 

27.236 


59 

2.298 

62.733 

2.901 

53.949 

1.911 

55.582 

2.086 

81.012 

13 

.014 

2T.653 

.011 

21-330 

.010 

16.057 

.014 

25.029 


60 

2.129 

63-578 

2. 61 3 

54.991 

1.778 

56.547 

1.938 

81 . 555 

14 

.020 

26.006 

.016 

19.906 

.015 

15.201 

.020 

23.163 











15 

.029 

24.350 

.024 

18.528 

.021 

12.34T 

.029 

21-720 


61 

2.007 

64.411 

2.651 

56.026 

1.653 

57.514 

1.798 

02.053 











62 

1.8T1 

65.23a 

2.495 

57.060 

1.535 

5B.491 

1.666 

62.517 

16 

. 04l 

22.689 

.034 

17.233 

.031 

9.625 

.042 

20.745 


63 

1.743 

66.064 

2.346 

58.099 

1.423 

59.479 

1.541 

82.954 

17 

.058 

21.027 

.048 

16.056 

.045 

7-307 

.060 

20 . 213 


64 

1.622 

66.893 

2.202 

59.146 

1.316 

60.482 

1.425 

83.370 

18 

.082 

19.356 

.067 

14.995 

.064 

5-883 

.084 

19.998 


65 

1.507 

67.730 

2.065 

60.206 

1.219 

61.503 

1.315 

83-770 

19 

-111 

17.678 

.092 

14.031 

.090 

5.667 

.146 

19.893 











20 

.148 

16.002 

.123 

13.137 

.122 

6.967 

.155 

19.666 


66 

1.399 

68.578 

1.934 

61.281 

1.127 

62.545 

l.a2 

84.160 











67 

1.297 

69.441 

1.810 

62.375 

1.0U0 

63.610 

1.115 

84.542 

21 

.192 

14.368 

.161 

12.286 

.162 

e.462 

.201 

19.145 


68 

1.201 

70.320 

1.691 

63.492 

• 959 

64.700 

1.027 

84.922 

22 

.241 

12.861 

.204 

11.467 

.208 

9.834 

.25 2 

18.278 


69 

1.111 

71.219 

1.578 

64.634 

.883 

65. 818 

.943 

85 . 301 

23 

■295 

11. 611 

.253 

10.691 

.260 

10.777 

.307 

IT. 169 


70 

1.027 

72.l4o 

1.472 

65.803 

.812 

66.966 

.866 

85.882 

24 

.351 

10.726 

.306 

10.001 

.316 

11.205 

.364 

16.121 











25 

.410 

10.222 

.362 

9.447 

.374 

11.190 

.425 

15-604 


71 

.948 

73.086 

1.370 

67.004 

.746 

68.146 

.793 

86.068 











72 

.874 

74.056 

1.275 

68.237 

.685 

69.360 

.726 

86.460 

26 

.475 

10.254 

.421 

9.104 

.438 

11.094 

.497 

17.80? 


73 

.804 

75.059 

1.184 

69.507 

.628 

70.609 

.664 

86.860 

27 

• 553 

11.765 

.486 

9.318 

-516 

11.681 

.592 

22.441 


74 

.740 

76.090 

1.099 

70.815 

• 575 

71.896 

.606 

87.269 

26 

.660 

15.552 

.574 

10.972 

.625 

13.756 

.726 

28.473 


75 

.679 

77.155 

1.018 

72.164 

• 525 

73.223 

• 553 

87.690 

29 

.612 

20.244 

.692 

14.380 

.790 

16.980 

.913 

32.753 











30 

1.028 

23.566 

.863 

18.175 

1.033 

19.715 

1.165 

33-801 


76 

.623 

78.254 

.942 

73.557 

.479 

74.592 

,504 

88.122 











77 

.571 

79-389 

.871 

74.995 

.437 

76.004 

.45B 

B8.567 

31 

1.319 

24.298 

1.103 

20.610 

1.565 

20.585 

1.479 

31.538 


78 

.522 

80.563 

.804 

76.483 

.396 

77.462 

.416 

89.027 

32 

1.684 

22.J47 

1.423 

20.637 

1.779 

19.339 

1.647 

27.862 


79 

.477 

81.778 

.741 

T6.023 

. 362 

78.969 

• 378 

89.501 

33 

2.110 

19-145 

1.810 

19-078 

2.248 

16.494 

2.250 

23-017 


80 

.435 

83.035 

.682 

79.617 

,328 

80.525 

,342 

89-992 

34 

2.574 

15.129 

2.275 

16.020 

2.736 

12.759 

2.668 

I6.4e2 











35 

3.046 

11.704 

2.767 

12.487 

3.207 

8.669 

3.080 

15.624 


61 

.396 

84.332 

.637 

81.269 

.298 

82.129 

.310 

90.498 











82 

.360 

85-665 

• 575 

82.983 

.270 

83.775 

.280 

91.010 

36 

3.499 

10.417 

3.267 

9.462 

3.630 

5.992 

3..4S7 

15.879 


63 

.328 

87.029 

.427 

84.744 

.244 

85-455 

.254 

91.554 

37 

3.906 

12.060 

3.747 

8.321 

3.964 

6.270 

3/815 

18.510 


64 

.298 

88.416 

.483 

86.537 

.2a 

87.163 

.230 

92.067 

38 ' 

4.256 

15.553 

4.182 

9-820 

4.259 

9.342 

4.116 

22.923 


85 

.272 

89-927 

.443 

86.358 

.200 

88.897 

.208 

92.611 

39 

4.535 

19.726 

4.557 

12.929 

4.454 

13.239 

4.363 

26.010 











40 

4.741 

24.047 

4.663 

16.594 

4.572 

17.259 

4.557 

33. 337 


86 

.247 

91.261 

.4o6 

90.207 

.181 

90.656 

.188 

93.164 











67 

.225 

92.718 

.373 

92.084 

.164 

92.442 

.170 

93.727 

41 

4.878 

28.309 

5.096 

20.354 

4.622 

21. 206 

4,6dO 

38.71? 


68 

.205 

94.199 

.342 

93.990 

.149 

94.253 

.154 

94.300 

42 

4.950 

32.433 

5.256 

24.068 

4.6l6 

25.018 

4.794 

44.0T4 


69 

.187 

95.702 

• 314 

95-924 

.135 

96.091 

.140 

94.883 

43 

4.967 

36.392 

5-356 

27.641 

4.562 

28.675 

4.847 

49- 38l 


90 

.170 

97.224 

.288 

97-887 

.123 

97.949 

.127 

95.473 

44 

4.937 

40.179 

5.396 

31.061 

4.473 

32.176 

4.663 

54.830 











45 

4,869 

43.798 

5.387 

34.323 

4.356 

35.528 

4.847 

59.824 
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distribution, and, consequently, the temperature distribution, would be 
changed. This approach could be a fruitful area for additional investi- 


gation when more data are available. 



V. COMPARISON OF THE MODEL WITH DATA 


The MRH data tape 

In order to estimate the various statistical parameters of the 
temperature distribution at each altitude as required in the temperature 
model, a computer tape containing experimentally measured temperature 
profiles v as obtained from the Aero-Astrodynamics Laboratory at the NASA 
Marshall Space Flight Center. The tape consists of kh2 sounding-rocket 
measurements of atmospheric temperature, pressure and density made at 
altitudes above 35 km by 25 Meteorological Rocket Network stations in 
different parts of the world. Of the 25 stations, seven are in the 30° 
latitude band, but only five of these had temperature profiles. (A 
number of stations reported only density or pressure profiles. ) These 
five stations and their longitudes are listed in Table 5*1* The columns 
under the heading "number of useable profiles" give the total number of 
independent daylight profiles at each site in the season indicated. In 
order to have even the small numbers listed here, "daylight" was defined 
as being between 5 a.m. and 7 p.m. If, for example, the winter daylight 
hours were redefined to be between 7 a.m. and 5 p.m. , six of the 20 
winter measurements would have been disqualified. 

In several instances two or more firings from a single station were 
made on the same day, sometimes only minutes apart. Since measurements 
cannot be considered independent in such circumstances, all temperature 
profiles made on the same day, at the same site, were first averaged to 
give a single profile which was then counted as a part of the sample of 
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Table 5.1 - Meteorological Rocket 



Latitude 

32.85° M 
30.38° N 
3^.12° N 
31.11° S 
2li.82° S 


Longitude 
106.10° N 
86.70° W 
119.12° W 
136.97° E 
113.87° E 


TOTALS 


Stations in 30° Latitude Band 
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independent profiles. 

Each profile consisted of a set of measurements made at one-kilo- 
meter intervals, but the range of altitudes covered varied from one 
profile to the next. Consequently, the number of observations at each 
altitude level varied as shown in Table 5*2. 

In addition to their use in estimating statistical parameters of 
the temperature distribution, the MRN data were used for comparison 
with model— generated temperatures, pressures and densities in an effort 
to eva l u a te the model. For each of the four seasons in the 30° latitude 
band, a sample of 500 pseudo-random atmospheres was generated (i.e., 

500 consistent sets of temperature, density and pressure profiles). 

The means and standard deviations of each property were computed at 
one-kilometer altitude intervals between 0 and 90 km. In addition, 
histograms of temperatures were constructed at various altitudes, and 
interlayer correlations among temperatures were computed. These model 
characteristics were then compared with corresponding characteristics 
of the MRN data. 

Comparison of model temperatures with data 

In order to compare trajectories computed using the present model 
with those based on the 1966 Standard Atmosphere Supplements, [2], it 
was decided to use temperatures prescribed by the latter as mean 
temperatures in the model. That is, the mean temperature at altitude 
z , for a particular season and latitude band, is defined as 



Table 5.2 - Number of MRN measurements at each altitude, 36 < z < 90 


Altitude 

Spring 

Summer 

Autumn 

Winter 

Altitude 

Spring 

Summer 

Autumn 

Winter 


0 

C 

0 

2 

64 

9 

5 

10 

15 


0 

0 

0 

2 

65 

9 

4 

9 

15 

38 

0 

0 

0 

3 






39 

1 

0 

1 

4 

66 

9 

4 

9 

14 

40 

4 

0 

3 

6 

67 

7 

3 

10 

l4 






68 

7 

A 

c 

8 

13 

4l 

4 

1 

4 

6 

69 

5 

1 

8 

12 

42 

6 

2 

6 

9 

70 

5 

0 

7 

11 

43 

8 

2 

9 

12 






44 

8 

2 

10 

13 ! 

71 

5 

0 

6 

11 

45 

8 

3 

10 

15 

72 

4 

0 

6 

12 






73 

3 

0 

6 

12 

46 

8 


13 

15 

74 

3 

0 

6 

12 

47 

8 

4 

13 

15 

75 

2 

0 

6 

12 

48 

9 

4 

13 

15 






49 

9 

4 

13 

17 

76 

2 

0 

6 

12 

50 

9 

4 

13 

17 

77 

2 

0 

6 

12 






78 

2 

0 

6 

12 

51 

9 

5 

13 

17 

79 

2 

0 

6 

12 

52 

9 

5 

13 

17 

80 

2 

0 

6 

11 

53 

9 

5 

13 

17 






54 

9 

5 

13 

18 

81 

2 

0 

6 

11 

55 

9 

5 

13 

17 

82 

2 

0 

6 

11 






83 

2 

0 

6 

10 

56 

9 

5 

13 

17 

84 

2 

0 

6 

9 

57 

9 

5 

13 

17 

85 

2 

0 

6 

9 

58 

9 

5 

13 

18 






59 

9 

5 

13 

18 

86 

2 

0 

6 

' 9 

60 

9 

5 

13 

18 

87 

2 

0 

6 

9 






88 

2 

0 

6 

9 

61 

9 

5 

12 

17 

89 

2 

0 

6 

7 

62 

9 

5 

12 

17 

90 

2 

0 

5 

7 

63 

9 

5 

11 

15 
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^(z) “ T^g (z; season, latitude) (5.1) 

where (zj season, latitude) is the temperature at altitude z 

oo 

given by the 1966 Standard Atmosphere Supplements, [2], for the season 
and latitude indicated. For the 30°, 60°, and 75° latitude bands, 
reference [2] provides only summer and winter models. In these cases, 
the spring/fall temperature profile is defined as the average of summer 
and winter profiles. 

Figure 5.1 compares model temperature means (solid curve) with 
the means of temperatures in the MRN data sample (circles). Ninety- 
five percent confidence intervals are drawn for each data point. The 
dashed curve on each plot is the temperature profile from the 1962 
Standard Atmosphere, [l], drawn here as a fixed reference to aid in 
comparing the four seasonal means. Seasonal variations are not as 
pronounced in the 30° latitude band as they are in more northern lati- 
tudes. 

The standard deviation of temperature at each altitude was estimated 
from the MRN data. As one might suspect from Figure 5-2, the amount of 
variation in the data estimates of a T (z) (circles) over the altitude 
range is within the amount expected for such small samples if one were 
sampling from populations with the same standard deviation. In other 
words, there is no significant difference between the standard deviations 
of temperatures at different altitudes in the MRN data. Therefore, it 
was assumed that, for a particular season-latitude model, the standard 
deviation of temperature is the same at all altitudes. Its value was 
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estimated by s 

P 

given by 


the usual pooled estimate of a standard deviation. 




(5.2) 


where n is the number of data points at altitude z ^ (see Table 5*2); 

s is the st an dard deviation of sample temperatures at altitude z 
l 1 

based on the MRN data, and the summation is only over i’s where n ± # 0. 

The four values of s used in the 30° latitude model (dashed lines on 

P 

Figure 5.2) are s (spring) - 7.45° K, s ( summer) - 6.82° K, s (autumn) - 
p P ■*' 

10,44° K, and s (winter) = 10.89° K. 

P 

Figures 5*3 - 5*6 show histograms of model-generated temperatures 
at lQ-km altitude intervals between 40 and 90 km. Each histogram 
represents a sample of 500 temperatures, and within each histogram 
class or subdivision, the circles represent data points lying in that 
interval. 

Perhaps the most important feature of the current temperature 
model - that which distinguishes it from any other model - is its inter- 
layer correlation structure (i.e., the correlations between temperatures 
at different altitudes). This correlation structure is linked directly 
to the behavior of the heating-rate profile T & (z) as affected by the 
ozone variation. Letting r^z^Zg) be the coefficient of correlation 
between T(* ) and T(Zg) , then in the present model 

r T (z 1 , z 2 ) = A 2 r s ( Z;L , z,,) 


(5-3) 











Relative frequency Relative frequency Relative frequency Relative frequency . Relative frequency 
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where r (z., z„) is the correlation between heating rates T (z_) 

Sl d. Si 

• * 

and T (z 0 ) , and X is the correlation between T(z) and T (z), 

3 cl S 

(see equation (3.7)). As atmospheric ozone varies according to its 

♦ 

distribution, the T profile follows a distinct pattern. The segment 

s 

of the T curve above its maximum at ~50 km altitude, (see Figure 

3 

1+.7), shifts to the right and left in such a manner that it produces 

4 

either a simultaneous increase or a simultaneous decrease in all T 

s 

values above 50 km. Thus, if and z^ are both greater than 50 

km, r s ^ z i» z 2^ approximately equal to +1, and hence ( z ^ , z^) is 

2 

approximately equal to X 

The solid curves in Figure 5.7 are plots based on the winter model, 
of r T ( z 1 * z 2 ^ versus z^ for nine values of between 35 and 90 

km. The winter model was chosen for comparison purposes because there 
are more MEN data for this season; plots of r,j(z^, * or °^^ er 

three seasons are quite similar. In the plots for * 60, 70, and 
80 km, note that r^Cz^, z^) , 50 < < 90, is approximately equal 

to X 2 « .81. 

» 

Below the stratopause, ~50 km, the behavior of T is more 

s 

* 

irregular. In the upper stratosphere (1+0 - 50 km), T values are 

S 

virtually uncorrelated to those above 50 km but have a strong negative 
correlation with heating rates in the lower stratosphere (20 - 30 km). 

This structure is reflected in the temperature correlations for z^ - 
1+0 and 1+5 km in Figure 5.7. 

The data symbols (circles) in Figure 5.7 are estimates of r T (z^, z^) 
based on the MRN data. The solid circles are those whose 95 percent 
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confidence intervale do not cover the model value of ^(z^, z 2^ * 

Because of the small sample sizes involved, it is difficult to make 
meaningful statements about the true nature of r^,(z^, 7.^) based on 
the data; nevertheless, there is reason for disappointment in the poor 
agreement seen in Figure 5*7* !h contradicting the model, these data 
estimates of r (z 1# z £ ) tend to imply that there is no strong linear 
correlation between T and T . One possible explanation for this is 

3 

that 10 of the 20 MKN profiles used to estimate r^z^ z^) were mea- 
sured either before 9 a.m. or after ^4 p.m. , times of day when the 
heating effects of the sun (particularly in winter) are not as pronounced 
as they are nearer to midday. Furthermore, the data represent varying 
sun zenith angles whereas the mo del- generated sample used a fixed 

zenith angle. Allowing this angle to vary in the model might produce 

* 

a different statistical distribution for the T g profile, thus 
affecting the distribution of temperatures. For these reasons, the 
comparison of model with data in Figure 5.7 is believed to be inconclu- 
sive. More data will be necessary in order to either validate or in- 
validate the model. 

The value of A = .9 used in this model determines the location 
of the model r^z^ z^) curves at points where z ± and z 2 are 
between 50 and 90 km. A smaller value of A , say A = .7 , gives 
better agreement with data points (circles) to the left of the curve 
where r^z^ z 2 ) - A 2 , but it also results in significant differences 
between model and data at virtually all points where |z 1 - z 2 | is 


small . 



Comparisons of model densities and pressures with data 


As explained in Chapter III, the natural logarithm of atmospheric 
pressure, log(P) , is the independent variable used as a vertical scale 
in the model. The reason for this is that the ozone model gives the 
partial pressure of ozone as a function of log(P) , and , hence, the 
solar heating rate profiles are based on this scale. 

To be more exact, the vertical scale used in the model is 

U = log(P/P Q ) 


where P is sea-level pressure. A set of fixed U values {U , 
o 

. . , U 9Q } is defined by 


u i = *«( p 66 ( ‘i )/P 66 (B o )) (5-5) 

where P^fz ) is the pressure at altitude z , (a. - 0, !***•» hm) , 
66 i 1 x 

as defined by the 1966 Standard Atmosphere Supplements, [2], for the 
appropriate season and latitude* Then at each a pseudo— random 

temperature is given hy 


[T (U ) - u (U. )] , „ 

T(U.) = y T (U. ) + Xa T ( u.) | (u -p + e <V (5*6) 

S 1 


where y^Clh) is the mean temperature at CL defined by 


V V = T 66 (a i ) 


(5*T) 


and 0^,(1^) “ 


is a constant estimated from the MRN data, as explained 
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above . 


From the T(U.) profile, one can then compute the geometric altitude 

0 


z(U. ) as 
1 


z( V 


1 , 

i 


T(u) du 


Mg 

~K 


- - f T(u) du 

° Ju. 


(5.6) 


where M, g Q , K, and R o are constants. (See Appendix A for derivation 

of 5.8). 

Next, using linear interpolation, the profile {T^), z(U i )} is 

converted to (T(z.), z } where the z.’s are even altitude increments 
11 

{0, 1, ...» 90 km} . 

Pressure and density profiles are then calculated for the even 
altitude increments as 


p( b ) = p b 



Mg. , 

KT dz 


(5.9) 


and 

M P(z. ) 

p ^ Z i^ = K tT^T ’ 

where P and z , in equation (5.9) are pressure and altitude, respect- 

b b 

ively, at some boundary point. It is customary to take - 0 and, 

hence P is sea-level pressure. However, since this is a statistical 
’ b 

model, one must specify the statistical distribution of P^ at the 
boundary altitude z^ , and then sample a value of P^ from that 
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distribution. It was found that, provided the same mean was used for 

P , the variance and shape of the P. distribution had negligible ef- 
b d 

fects on the means of pressures and densities. These choices did, how- 
ever, have a fairly sizable impact on the variances of pressure and 
density. The selection of and its associated P^ distribution will 

be discussed in more detail shortly. 

Figure 5.8 shows mean density profiles (solid curves) corresponding 
to the four 30° latitude models. The actual quantity plotted here is 
the percent departure from the 1962 Standard density P^g^ 2 ) as ^ e -^ net ^ 
by 


V 2) 


P(z) - P 62 (z) 
P 62 U) 


x 100 % 


(5*H) 


where p(z) is the average of 500 pseudo-random density values generated 
by the model at altitude z . Since densities and pressures vary by 
about six orders of magnitude between sea level and 90 km, pressures 
and densities are customarily plotted as a percent departure from the 
1962 Standard pressures and densities, thus avoiding the use of 6— cycle 
log paper. 

The dashed lines on Figure 5*6 represent the 1966 Standard densities 
for the four seasons, the density profiles which correspond to the mean 
temperature profiles used in the model. As one can see, there is close 
agreement between the 196 6 Standard and the model's mean densities. 
However, perfect agreement between the Standard density and the model's 
mean density should not be expected. Indeed, such agreement would be 
surprising, since, in general, if F is a function of a random variable 




Altitude . km 



T9 


X , the mean of which is P x , the mean of F(X) is not necessarily 
equal to F(y^). 

The data symbols (circles) on Figure 5.8 show the percent departures 
of the MRN density averages from the 1962 Standard density. Again 95 
percent confidence intervals are drawn. At best one can say that the 
data averages and model averages exhibit the same seasonal trends. 

Figure 5.9, similar to Figure 5*8, is a comparison of the percent de- 
partures of pressure means from the 1962 Standard Pressure. 

As stated earlier, the choice of the distribution at the 

boundary point has a significant impact on the variances of pres- 

sures and densities. Consequently, this choice could not be arbitrary, 
but had to be justifiable in terms of physical reality. Accordingly, 
it was decided to select as a boundary point an altitude where density 
variation is minimal. There are certain altitudes, called isopynic 
layers, [1*0], where atmospheric density is almost constant throughout 
the year. These isopynic layers lie at altitudes of approximately 8, 

2U, and 90 km. By taking z^ to be one of these altitudes and by 
letting p.^ , atmospheric density at altitude , be a constant, 

then the distribution of P is determined by the distribution of 


T(z ) since the equation of state gives 
b 


p t = 


K T<z t ) 


( 5 . 12 ) 


Two boundary points , = 2l km and z^ = 90 km , were used in 

the present model. When z. < 80 km, in equation (5.9), the boundary 
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point - 24 km is used, and when z. > 80 km, ^ = 90 km is used. 

The selection of 80 km as a dividing point was made because at z ± - 80 
km the variance of atmospheric density is the same using either 2.H 

or z. = 90 km. 

D 

Although it was only necessary to select one boundary point , the two 
were chosen here to force the existence of isopynic layers at those 
altitudes. The 8 km isopynic layer, on the other hand, has been omitted 
since it does not affect the 30 - 90 km region of primary interest. If 
the model were perfect, then isopynic layers would be created naturally 
at 8, 24, and 90 km even if only one of these is used as the boundary 
point. However, in the present model, the use of only one boundary 
point resulted in the disappearance of isopynic layers at non-boundary 
points . 

In Figure 5. 10 the ratio of the standard deviation of density to 
the 1962 Standard density is plotted for each season. The solid curves 
represent standard deviations calculated from samples of 500 model- 
generated density profiles. The data symbols are the corresponding 
standard deviations in the MRN samples, and 95 percent confidence 
intervals are drawn. Note the isopynic layers at z = 24 km and 90 
km. The peak in standard deviations at 80 km is a result of using this 
altitude as a dividing point, as explained above, but the peak at around 
30 km results solely from the way the temperatures vary and is not re- 
lated to the integration boundaries. 

Figure 5.11 is analogous to Figure 5.10 except that here standard 
deviations of pressure are plotted. The solid curves, representing the 







.2 .3 .4 
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deviations of pressure with 


of pressure. 
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model standard deviations, resemble those on Figure 5*10. However, the 
minimum points appear to be approximately 5 It® below the isopynic levels, 
and the peak at approximately 35 km is less pronounced. There is a 
discontinuity in the standard deviations of model pressures at 80 km, 
which is caused by the use of 80 km as an integration boundary. 

Conclusions 

In comparing the model with data, a distinction should be drawn 
between model characteristics which are estimated from the data and 
those which are not. Temperature means and standard deviations are 
examples of the former, whereaB , the statistical properties of pressure 
and density are of the latter type. 

Since temperature means and standard deviations are selected by the 
model user, the amount of disagreement between these and their counter- 
parts in the data is a matter of choice. It was decided, for example, 
to use the 1966 Standard Atmosphere temperatures as model means instead 
of estimating means from the data. In the case of temperature standard 
deviations, a constant standard deviation was used for all altitudes, 
and this value was obtained by pooling standard deviations in the MEN 
data sample. Thus, any differences between model and data seen in 
Figures 5.1 and 5.2 do not constitute errors in the model. 

The means of model pressures and densities agree well with their 
corresponding 1966 Standard Atmosphere values even though such agree- 
ment was not guaranteed a priori. As in the case of temperature, the 
data estimates of pressure and density means differ from the 1966 



Standard and, therefore, cannot be expected to match model means. 

The standard deviations of model pressures and densities, which 
ideally should agree with data about as well as do temperature standard 
deviations, instead show slightly poorer agreement. If the integration 
scheme could be modified so as to remove the artificial peak at 80 km, 
agreement between model and data standard deviations would be signifi- 
cantly improved. 

The most serious discrepancy between model and data lies in the 

I 

interlayer temperature correlations, r^(z^, z^) . Although some of 

this error can perhaps be explained in terms of zenith-angle differences, 

it still may be necessary to alter the model. One possible solution 

would be to extract from the unmodeled error term a source of variation 

which correlates temperatures at adjacent altitudes. Winds would be 

such a source, for example. Thus, X could be reduced and most model 

r ( 2 z ) values made smaller without decreasing the correlation be- 
T 1 2 

tween temperatures at adjacent altitudes. 

Another possible method of improving the agreement between the data 
and model r (z^, z^) values would be to allow X to vary as a function 
of altitude. There is physical Justification for this since X , the 
correlation between temperature and the solar heating rate of ozone, 
is smaller at altitudes where there is absorption by other gases (e.g. , 
water vapor) than where 0^ is the only absorber. Before attempting 
to improve the agreement between model and data, however, the data 
base should be enlarged so as to increase the confidence in estimates 



VI. THE EFFECT OF ATMOSPHERIC VARIABILITY 
OH FIVE IMPORTANT SPACE SHUTTLE ENTRY PARAMETERS: 

AH APPLICATION OF THE MODEL 

A Monte Carlo entry trajectory program was developed "by modifying 
a rapid, three-dimensional point— mass entry program to use pseudo-random 
atmospheric density profiles instead of the 1962 U. S. Standard density 
profile, [l]. Large samples of entry trajectories were then generated, 
and the statistical behavior of five major performance indicators was 
studied. 

The five entry parameters selected for study were: 

(1) maximum dynamic pressure 

(2) maximum stagnation— point heating rate 

( 3 ) maximum g-load 

(ii) final down-range distance 

(5) final cross-range distance 

The objective of the study was to analyze the statistical distributions 
of these five parameters and, thereby, determine "three-sigma” design 
values for each. The adjective "three-sigma" is used here to refer to 
a value which will be exceeded with a probability of .0013. If the 
statistical distribution of the parameter is normal (Gaussian), then 
this design value will, in fact, lie a distance of three standard de- 
viations (three sigma’s) away from the mean. In general, however, 
this .0013-probability design value is not three standard deviations from 
the mean, and hence, the quotation marks on "three-sigma". 
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A two-dimensional random density model 

Instead of entering along a steep trajectory easily described by 
one atmospheric profile, the Space Shuttle will actually fly back from 
orbit along a shallow trajectory which takes it a distance of several 
thousand kilometers down range. (See Figure 6.1.) Thus, a stochastic 
atmosphere model used to simulate Shuttle entries should provide vari- 
ation in atmospheric properties along horizontal flight paths i one 
random atmospheric profile per entry assigning a single set of properties 

at each altitude is not adequate. 

In the present simulation the model described in the preceeding 
chapters was used to generate one pseudo— random atmospheric density 
profile at the location where entry began. Density profiles at locations 
down-range differ from this initial profile but are not completely un- 
correlated. Therefore, at a location 110 km down range from the initial 
entry point, it would have been wrong to generate a new independent 
atmospheric density profile using the same model. Instead, the initial 
profile was perturbed slightly to give a new profile 110 km down range. 

An explanation will be given below as to how this was done and why the 
110 km increment was used. At a point 220 km down range a new density 
profile was obtained by perturbing the previous profile, (i.e., the one 
at 110 km down range), and this was continued until a sequence of density 
profiles spaced 110 km apart was generated. The Shuttle orbiter was 
then flown through this array of densities, and the value of atmospheric 
density at any point along its trajectory was found by linear interpolation 




Down range 



Figure 6.1 - Illustration of shallow entry path of Space Shuttle Orbiter. 
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"between nearby density profiles. 

Let p (z) denote a pseudo-random value of atmospheric density at 

o 

altitude z (z = 0, 1,...» 90 km) generated by the stochastic atmosphere 
model described in the preceeding chapters. The profile {p Q (z): z 

0, 90 km} is located at the point where down-range distance DR 

is 0. Let { p ( z ) } be the density profile located where DR = 110 km, 
and in general, let {p^ (z )} be located where DR = J x no km. Then 
the density Pj + -^( z ) was gi ven t>y 

p (z) = + ^ z ) J = 1 » 2 » •** (6.1) 

j+l J 

where 

6(z) - N(0, (z)) (6.2) 

The variance of 6(z) , a^(z) , was based on a set of maximum hori- 

zontal density gradients given in reference [23]. This reference 
places an upper limit on the horizontal change in atmospheric density 
between two locations 110 km apart (hence the use of 110 km intervals). 

Since p (z) and p (z) » as defined above, are two densities at 
j J ^ 

altitude z located 110 km apart, then according to [23] the following 
must hold: 

-y(z)p 62 (z) < p j+ 1 u) - pj( z > i yU)p 62 ( z ) (6 - 3) 

where Pg 2 ( z ) is the 1962 Standard atmospheric density at altitude z , 
and y( z ) is obtained by linear interpolation from Table 6.1. 

In the present situation p^ +1 (z) - p^(z) = M z ) is a normal 
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Table 6.1 - Design Horizontal Density Gradients 


30° Latitude Band 


Altitude 

Spring/ Autumn 

Summer 

Winter 

0 

0.35 

0.20 

0.50 

30 

0.35 

0.20 

0.50 

hO 

0.70 

0.30 

1.10 

50 

0.95 

0.40 

1.50 

60 

1.25 

0.50 

2.00 

TO 

l.Uo 

0.60 

2.20 

80 

1.25 

0.50 

2.00 

90 

0.20 

0.10 

0.30 

If 

p and p p are two values of 

J- Cm 

density at the same 

altitude z 

but 110 

km apart, then |p^ - p^l must 

be less than TPg 2 

, where y 


is interpolated from this table, and is the 1962 Standard Atmo- 

sphere density at altitude z . 


Table 6.2 - Shuttle Characteristics used in Simulations 


Orb iter weight 

= 

870000 N (196000 lbs) 

Reference area 

= 

268 m 2 (2890 ft 2 ) 

Angle of attack 

= 

3*+ deg 

Lift-to-drag ratio 


1.25 

Drag coefficient 

= 

0.572 

Lift coefficient 

= 

0M1 


2 ° 


random variable and, hence, does not have strict upper and lower limits 

as required by reference [23], equation (6.3). However, by chosing 

o.(z) such that 
o 

Y(z)P 62 (z) = 3o a (z) < 6 - U 

the absolute limit in equation (6.3) was replaced by 

Pr [-y(z)P^ 2 (z) < Pj +1 ( z ) “ P J^ - Y ^ P 62^ Z ^ = 

a probabilistic limit. The design horizontal density gradients of 
re f erence [23] then became "three- sigma” design values. 

Procedure 

The vehicle and aerodynamic characteristics used in this study 
(see Table 6.2) are based on the current Shuttle design, [4l]. Four 
samples of 1000 pseudo-random entries were generated, each sample 
representing a different season in the 30° latitude band. Within 
each sample, the only difference between any two entry trajectories 
was their respective atmospheres. The five entry parameters listed 
above were recorded for each entry. Sample means and standard devia- 
tions were calculated; histograms were constructed, and empirical 
cumulative frequencies were plotted on normal probability paper. The 
first four moments of each parameter were used to fit a Pearson 
distribution, [*+2], to the empirical distributions (see Appendix D). 
For comparison purposes, a number of different methods were used 
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to determine "three-sigma" design values. The methods used were the 
following : 

Method A (Traditional Method) . The parameter was assumed to have a 
normal distribution, and its design value was given by 

(A) design value = X t 3s (6.6) 

where X and s are the sample mean and sample standard deviation, 
respectively. A plus or minus sign was chosen depending on the nature 
of the parameter. 

Method B (Probability Paper Method) . Again, the parameter was assumed 
to have a normal distribution. However, by this method, a normal 
cumulative frequency curve F(x) was fitted to the empirical cumulative 
frequencies, using probability paper, and the design value was the 
parameter value, x = (B) design value, satisfying F(x) = .9987 
(or .0013 if it is a minimum). 

Method C (Non-parametric Method) . The "three-sigma" design value , 
as defined above, is a value which will be exceeded approximately once 
in a sample of size 1000. Letting x^ and x^ ^ denote the largest 
and next largest values, respectively, in the parameter sample of size 
1000. then any value between x _ and x has been exceeded once in 
that sample and could justifiably be chosen as the "three-sigma" design 
value. Method C took the design value to be 

(C) design value = x n 


( 6 . 7 ) 



92 


the most conservative choice between x n * ^ ln Z 

where a design value was a minimum the smallest sample value was used). 

Method D (Pearson Distribution Method) . By this method, the first four 
sample moments were used to fit one of the Pearson distributions to the 
observed distribution. Then the "three- sigma" design value was that 
value which satisfied F p (x) = .9987 (or .0013) where F p (0 was the 
cumulative distribution function of the Pearson distribution which best 
fitted the data. 


Results and Discussion 


Table 6.3 lists various sample statistics computed from the four 
entry samples. Values in the "Deterministic Prediction" column are 
values of the parameters obtained from non-random trajectories using 
the appropriate atmosphere from the 1966 Standard Atmosphere Supplements, 

[ 2 ]. 

The coefficients of variation, defined by. 


CV 



x 100# 


( 6 . 8 ) 


can be used to estimate confidence intervals for the means. A 95 percent 
confidence interval for any mean in Table 6.3 is 

X (1 ± .0006 CV) (6.9) 


For example, since the coefficient of variation of autumn down-range 
distances is 3#, one can state with 95 percent confidence, that the 
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Table 6.3 - Entry Parameter Statistics Based on 1000 Entries 


Parameter 


Deterministic Mean Standard 

Prediction Deviation 


Coefficient of 
Variation, % 


(a) Maximum dynamic pressure, 


N/m 2 


Spring 

Summer 

Autumn 

Winter 


5305 

5505 

5296 

5 U 31 

5305 

5592 

5307 

567 U 


104.1 

1.9 

86.4 

1.6 

158.3 

2.8 

168.5 

2.9 


(t) Maximum stagnation— point heating rate, 


x 10~ 5 , 


W/m £ 


Spring 

Summer 

Autumn 

Winter 


9.910 

10.059 

9-910 

9 . 71+9 


9.897 

10.069 

9.915 

9.863 


0.353 

0.375 

0.462 

0.4l4 


3.6 

3.7 

4.7 
4.2 


(c) Maximum 

Spring 

Summer 

Autumn 

Winter 


g-load 

1.165 

1.150 

1.165 

1.152 


1.206 

0.0251 

2.1 

1.185 

0.0180 

1.5 

1.223 

0.0346 

2.8 

1.247 

0.0412 

3.3 


(d) Down range, km 

Spring 13610 
Summer 13450 
Autumn 13610 
Winter 13790 


13630 

294.1 

2.2 

13470 

396.3 

2.9 

13650 

415.2 

3.0 

13760 

• 394.1 

2.9 


(e) Cross range. 

km 

Spring 

1290 

Summer 

1295 

Autumn 

1290 

Winter 

1283 


1288 

22.46 

1.7 

1297 

26.80 

2.1 

1286 

32.87 

2.6 

1279 

31.56 

2.5 



estimated mean down-range distance in autumn, X = 13650 km, is accurate 

to within i .0018 X (i.e., to within ± .18$). 

Figures 6.2 - 6.6 show histograms of the five entry parameters. 

The dashed lines are the Pearson probability density functions fitted 
to each distribution. One. general characteristic of each parameter is 
that its autumn and winter distributions have more dispersion (greater 
variance) than its spring and summer distributions. This reflects the 
fact that the variances of temperature, density and pressure in the 
autumn and winter atmosphere models were greater than those of the spring 
and summer models. The down-range distributions appear to be skewed 
to the left, whereas all other distributions appear to be skewed to the 

right . 

In the cases of the summer maximum dynamic pressure, down-range, 
and cross-range distributions, it was not possible to find Pearson 
distributions using the usual moments fit. This is because 3-^ 

0 were inadmissible (see definitions in Appendix D). It was possible 
to modify the B and $ 2 values for maximum dynamic pressure to get 
the "closest" Pearson fit. (See Appendix D). This is the dashed curve 
in Figure 6.2 (summer). The same modification when applied to the 
summer down-range and cross-range distributions resulted in J-shaped 
Pearson Type VI distributions which did not resemble the sample 
distributions at all. Therefore, no Pearson fits were obtained for the 
summer down-range and cross-range distributions. *' 

Table 6.1+ lists the various "three-sigma" design values obtained 
by Methods A, B, C, and D, as described above. If a design value lies 










Relative frequency . 
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Table 6.4 - Parameter Design Values 


A B C D 


Parameter 

X±3s 

Probability 

Paper 

Sample 

Extreme 

Pearson 

Distribution 

(a) Maximum dynamic pressure, 

N/m 2 



Spring 

5817 

5903 

6192 

5997 

Summer 

5690 

5772 

6131 

5913 

Autumn 

6067 

6201 

6298 

6310 

Winter 

6180 

6305 

6421 

6345 

(b) Maximum stagnation-point 

—6 

heating rate, x 10 , 

W/m 2 


Spring 

1.096 

1.126 

1.226 

1.161 

Summer 

1.119 

1.150 

1.192 

1.174 

Autumn 

1.130 

1.169 

1.204 

1.192 

Winter 

1.111 

1.150 

1.317 

1.225 

(c) Maximum g-load 





Spring 

1.281 

1.302 

1.325 

1.312 

Summer 

1.239 

1.254 

1.273 

• 1.265 

Autumn 

1.327 

1.355 

1.390 

1.377 

Wint er 

1.371 

1.401 

1.421 

1.409 

(d) Maximum down range 

, km 




Spring 

14510 

14680 

14310 

14260 

Summer 

14650 

14830 

21830 

# 

Autumn 

14900 

15090 

14780 

14610 

Winter 

14940 

15230 

14800 

14690 

(e) Minimum down range 

, km 




Spring 

12750 

12550 

11900 

12420 

Summer 

12280 

12200 

12140 

* 

Autumn 

12410 

12200 

12060 

12130 

Winter 

12580 

12330 

11450 

12100 

(f) Maximum cross range, km 




Spring 

1355 

137 1 * 

1431 

1393 

Summer 

1377 

1398 

1407 


Autumn 

1385 

1407 

l4l6 

1415 

Winter 

1374 

1400 

1482 

1428 

(g) Minimum cross range, km 




Spring 

1221 

1223 

1233 

1230 

Summer 

1217 

1224 

905 

« 

Autumn 

1187 

1191 

1203 

1205 

Winter 

1184 

1191 

1206 

1208 
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in the direction of skewness (i.e., in the long tail), then the tradi- 
tional use of X ± 3s as a design value will, in general, lead to an 
exceedance probability (risk) higher than the intended .0013 . Because 
of the positive skewness in the maximum dynamic pressure, maximum 
stagnation-point heating rate, and maximum g-load distributions, it is 
not surprising, therefore, that the (A) design values (X ± 3s) are 
less conservative (more "risky” ) than the (B) , (C) , or (D) design 
values which are based on sample cumulative frequencies. In other 
words, the latter three methods estimate the true cum ul ative distribution 
function F(x) and select a design value x satisfying F(x) = .9987* 
whereas the (A) design value, x^ , generally has a cum u lative frequency 
F(x, ) less than .9987* Taking the autumn maximum stagnation-point 
heating rate distribution as an example, and assuming the Pearson 
distribution shown in Figure 6.3 is the true distribution, the probability 
of exceeding the (A) design value (11.301), is .0119 instead of .0013 
as desired. The chance of exceeding this design value at least once 
in 100 missions (the anticipated Shuttle lifetime) is 70 percent, 
whereas the chance of exceeding the (D) design value (11.919) at least 
once is only 13 percent. 

In the cases of the minimum cross-range and the maximum down-range 
design values, the traditional Method (A) is actually more conservative 
than the (D) method. This is because these design values lie in the 
short tails of their respective skewed distributions. Nevertheless, 
over-conservatism is still not a desirable trait since it adds to the 


cost of the design. 
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If it is decided that the optimum exceedance probability for a 
design value should be .0013, and both greater and smaller exceedance 
probabilities are to be avoided, then the use of a Pearson fit is 
recommended whenever possible. Accordingly, the (D) design values 
listed in Table 6.3 are preferred to those given by the other methods. 

A pilot study similar to the present study was made by the author , 
[1+3], to determine if random fluctuations in atmospheric density, on the 
order of those observed, could have an appreciable effect on Shuttle 
entry parameters. (The same five parameters were studied). In that 
study a fairly crude stochastic atmosphere model was used. However, 
in spite of its deficiencies, the pilot study did produce some useful 
information. It showed, for example, that the shape of a parameter's 
distribution persisted over a wide range in entry angles of attack, 
and the coefficient of variation for a particular parameter remained 
quite constant even though the magnitude of its mean varied considerably 
for different angles of attack. 

There is every reason to believe that these same consistencies 
would apply in the present study. Therefore, even though angles of 
attack or other trajectory characteristics may change, thus changing 
magnitudes of the entry parameters, the shape of a parameter's distri- 
bution and its percentage variation due to atmospheric effects should 
not be seriously affected. For this reason, design values expressed 
as ’’design factors” to be multiplied by the parameter's mean are 
perhaps more valuable than strict design values. Table 6.5 lists 
these factors based on the (D) design values in Table 6.4. The latter 
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Table 6.5 - Parameter Design Factors BaBed on Pearson Distributions 


Parameter 


Factor 


(a) Maximum dynamic pressure 

Spring 

Summer 

Autumn 

Winter 


1.09 

1.09 

1.13 

1.12 


(b) Maximum stagnation-point heating rate 

Spring 

Summer 

Autumn 

Winter 

(c) Maximum g-load 

Spring 

Summer 

Autumn 

Winter 

(d) Maximum down range 

Spring 

Summer 

Autumn 

Winter 

(e) Minimum down range 

Spring 

Summer 

Autumn 

Winter 

(f) Maximum cross range 

Spring 

Summer 

Autumn 

Winter 

(g) Minimum cross range 

Spring 

Summer 

Autumn 

Winter 


1.17 

1.17 

1.20 

1.2b 


1.09 

1.07 

1.13 

1.13 


1.05 

1.09 

1.06 


0.91 

0.89 

0.88 


1.08 

1.10 

1.12 


0.95 


0 . 9 ^ 

0.9b 



are obtained by multiplying the design factors in Table 6.5 by the 
parameter’s mean given in Table 6.3. 

Conclusions 

This study has shown that significant variations can occur in entry 
parameters as a result of day-to-day atmospheric variability. Parameter 
distributions showed seasonal trends in both their means and th&tr vari- 
ances. Largest and smallest means occurred in the summer and winter 
distributions and intermediate means in the spring and fall. Autumn 
and winter variances were larger than spring and summer variances in all 
parameter distributions because the autumn and winter atmosphere models 
were more variable than the spring and summer models. 

Using the design factors in Table 6.5, general conclusions can be 
drawn about the relative a mo unt of variation one might expect to result 
from atmospheric variability. Design values for maximum dynamic pressure 
and maximum g-load range from 7 to 13 percent higher than nominal (mean) 
values; design values for maximum stagnation-point heating rate range 
between 17 and 2k percent higher than their nominal values; down-range 
design values are between +9 and -12 percent, and cross-range design 
values are between +12 and -6 percent of their nominal. 

Interpretation of the design values for maximum dynamic pressure 
and maximum g-load is straight-forward. The significance of maximum 
stagnation-point heating rate, however, is its impact on the total 
integrated heat load and its impact on the maximum stagnation-point 
surface temperature, T . The total heat load was not recorded 
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during the simulations, and hence, its variation is not known. A 
surface temperature at any point on the vehicle is proportional to 
the fourth root of the heating rate at that point. Therefore, 
is proportional to the fourth root of the maximum stagnation-point 
heating rate. This means that a design value for the latter which is 
24 percent higher than its nominal is equivalent to a design maximum 
temperature approximately 6 percent higher than its nominal (i.e., 1.24 
= 1.06) . ■ 

Departures from the nominal down-range and cross-range distances 
are not to he interpreted as miss distances since on-hoard guidance 
systems will allow the Shuttle orhiter to correct for flight path 
errors. The down-range and cross-range design values obtained here 
imply that these guidance systems must he capable of correcting for 
errors of as much as 12 percent caused by natural variations in at- 
mospheric density. 

The present study has shown that atmospheric variability is an 
important source of error to consider in Shuttle entry design studies. 

It is not the only error source, on the other hand, since uncertainties 
in the vehicle* s weight, aerodynamic coefficients and initial orbital 
elements will all cause errors in the entry trajectory. Past error 
analyses have often either ignored atmospheric variability or else 
used '* extreme" atmospheres as described in Chapter I. It is recommended 
that future error analyses include atmospheric variations and that a 
statistical approach such as the present one be used. 



VII. SUMMARY AMD CONCLUSIONS 


Summary 

A stochastic model of the atmosphere between 30 and 90 km was 
developed for use in Monte Carlo Shuttle entry studies. The model is 
actually a family of models , one for each latitude- season category as 
defined in the 1966 U, S. Standard Atmosphere Supplements, f 2 ] , 

Each latitude- season model generates a pseudo-random temperature 
profile whose mean is the appropriate temperature profile from the 1966 
Standard Atmosphere Supplements. The standard deviation of temperature 
for a given latitude-season model is assumed to be the same at all 
altitudes and is estimated from sounding-rocket data. Departures from 
the mean temperature at each altitude are produced by assuming a linear 
regression of temperature T(z) on T^(z) , the solar heating rate of 
ozone. First, a profile of random ozone concentrations is generated 
using an auxiliary stochastic ozone model, which was dlso developed as 
a part of this study. The solar heating rate T g is random since it 
is a function of the random ozone concentrations. The steps taken in 
generating a random temperature profile are illustrated schematically 
in Figure 7*1- 

Pressure and density profiles are calculated from each temperature 
profile by solving simultaneously the hydrostatic equation and the e- 
quation of state. Thus, each pseudo-random atmosphere consists of a 
consistent set of temperature, pressure and density profiles. Isopynic 
layers, regions where atmospheric density is almost constant, are created 
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Altitude . km 
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STEP \. Generate a random 
0 3 profile (solid line) 


STEP 2 . Calculate the solar heatinq 
rate (solid line) corresponding to 
random 0^ profile 



STEP 3. Adjust temperature protile accordinq STEP 4. Add secondary random fluctuation 
to whether solar heatinq rate is above or ' n0,se ' 

below its mean 


Figure 7.1 - Steps taken in generating one pseudo-random temperature profile 
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at 24 and 90 km altitudes by using these two altitudes as boundary points 

in integrating the hydrostatic equation. 

A computer data tape of 442 sounding-rocket measurements of the 
atmosphere above 35 km was obtained to use in estimating parameters 
of the temperature distribution. These data were also used for compar 
ing their statistical characteristics with those of the model. After 
the data were sorted according to latitude- seas on categories and adjust- 
ments were made to assure independence of profiles, the sizes of resul- 
tant samples were quite small. For example, in the 30 latitude band, 
sample sizes ranged from 5 id the summer sample to 20 in the winter 
sample. Because of such small samples, confidence intervals associated 
with the sample statistics (e.g., means, standard deviations, and 
correlations) were so wide that comparison of these statistics with 

model statistics were inconclusive. 

The four 30° latitude models were used in Monte Carlo simulations 
of Shuttle entries to study the effect of atmospheric variability on 
five Shuttle entry parameters. The simulations were of point-mass 
trajectories, and the five parameters studied were maximum dynamic 
pressure, maximum stagnation-point heating rate, maximum g-load, final 
down-range distance and final cross-range distance. A sample of 1000 
entries was generated for each of the four seasons, and the resultant 
parameter distributions were analyzed to determine design values. Pearson 
distributions were fitted to the empirical distributions, and design 
values obtained from these distributions were compared to the traditional 
’’three- sigma" design values which are based on assumptions that parameter 
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distributions are normal. Pearson design values were defined as those 
parameter values with exceedance probabilities equal to .0013 based 
on the fitted Pearson distribution. (The value .0013 is the exceedance 
probability associated with a normal random variable located a distance 
of three standard deviations from its mean. ) 

In general, the design values based on fitted Pearson distributions 
were more conservative (farther from the mean) than those obtained 
with the assumption of normality. This is because all of the empirical 
distributions were skewed, and most design values lay in the direction 
of skewness (long tail). Pearson design values ranged from being 2k 
percent greater than the mean, in the case of winter maximum stagnation- 
point heating rates, to being within 5 percent of the mean, in the case 
of spring minimum cross-range distances. 

Autumn and winter parameter distributions showed more dispersion 
than spring and summer distributions, because the autumn and winter 
atmosphere models were more variable than the spring and summer models. 

It is apparent from this fact that parameter design values are quite 
sensitive to the estimates of variances used in the atmosphere models. 

In particular, two factors contribute to the observed dispersion dif- 
ferences between summer and winter histograms. One is the difference 
between o T ( summer) = 6.82° K and (^(winter) = 10.89° K, and the other 
is the difference between summer and winter design density gradients 
(Table 6.1). Because of this sensitivity. Shuttle parameter design values 
are only as reliable as these estimates. This further underscores the 
need to increase the data base used in estimates of atmosphere model 
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parameters. 

Concluding observations 

There are two basic philosophies one can use in developing a 
stochastic model of some random phenomenon. By far the most common 
approach is to observe the available data relating to that phenomenon 
and then model that data empirically without any regard to cause-and- 
effect relationships. By using the data to estimate means, variances, 
correlations, and higher moments, one can generate artificial data 
samples with all the statistical characteristics of the original sample. 

A second alternate approach is to attempt to identify and model 
the fundamental random processes which bring about the observed random- 
ness in the phenomenon being modeled. For example, in the case of 
atmospheric temperature, one would look for underlying random processes 
which cause temperature to vary. 

Believing that a good theoretical model is preferable to a good 
empirical model, if both were attainable, the author chose the theoret- 
ical approach in setting out to model atmospheric variations. Although 
the development of a workable theoretical model is a much greater 
challenge, such a model would provide physical insight into the nature 
of atmospheric variations. A purely empirical model, on the other hand, 
is strictly limited to describing a particular set of data, and its 
user must be cautious when projecting the model beyond the range of 
the original data. Furthermore, empirical models generally are time 
consuming and require the storage of large correlation matrices. 



Ill 

The fact that the microscopic nature of any gas is basically 
random was further encouragement for taking a theoretical approach 
to modeling the atmosphere. For example, the temperature of a gas 
is actu ally a direct measure of the variance of molecular velocities 
in that gas. That is, the velocity of a molecule of gas in any one 
direction is a normal random variable with mean 0 and variance equal 
to KT/M, where K is the universal gas constant, M is the molecule's 
weight, and T is the gas temperature, [44]. 

In accordance with the "theoretical* 1 philosophy , one fundamental 

random process, the creation and destruction of atmospheric ozone, was 

identified as the major cause of temperature variation in the 30-90 km 

region* An attempt to model it on a microscopic level, using collision 

rates, etc., lead to the immediate realization that variations in a 

macroscopic gas property could not be accounted for by modeling random 

events occurring on a microscopic scale. The law of large numbers 

steps into play so that, if N is the number of molecules in the gas 

volume, the collective effect of N random events is virtually non- 

20 

random at altitudes below 90 km where N is greater than 10 molecules 
per cubic meter. 

Thus, the initial effort to model the causes of temperature variation 
using purely theoretical considerations failed. It led, instead, to an 
empirical model of atmospheric ozone based on experimental measurements 
of ozone concentrations, [37]* Pseudo-random ozone profiles were used 
to calculate the solar heating rate of ozone at each altitude, and these 
in turn were linked to temperatures via a linear regression model. 
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A reasonable question at this point might be, "Why avoid the use 
of an empirical temperature model and at the same time resort to the 
use of an empirical ozone model?" Admittedly, this is a compromise, 
and a theoretical model for ozone variation would be preferable. 

However, even though the present model contains no insight into the 
nature of ozone variations, it still uses these to explain, at least 
partially, the variations in temperature, and particularly the corre- 
lations between temperatures at different altitudes. The particular 
ozone model used here determines the shape of the temperature distribu- 
tion at each altitude. However, any ozone model showing the same uniform 
increase and decrease in ozone above 30 km, as exhibited by the data, 
will give the same interlayer correlation structure for temperatures ; 
it is this feature - the correlation structure - which largely governs 
the pressure and density distributions. Thus, the pressure and density 
distributions and the interlayer temperature correlations are somewhat 
robust to the choice of ozone models. 

A number of improvements can and should be made in order to obtain 
better agreement between the temperature model and data. A better, 
more quantitative, means of estimating X should be devised; other 
sources of variation should be included (e.g, , winds, water vapor ab- 
sorption, etc.), and a more acceptable meanB of estimating boundary 
conditions in the hydrostatic equation should be sought. A prerequisite 
to any improvement attempt, however, is the enlargement of the data 
base. As it stands now, the model represents an encouraging "first cut" 
at a workable theoretical model. 
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APPENDIX A. EQUATIONS FOR CALCULATING ATMOSPHERIC PROPERTIES 


List of Symbols 


Symbol 

8 

g o 

K 

M 

P 


Definition 


acceleration due to gravity 


Units/Constants 

m/sec 2 


acceleration due to gravity at sea level m/sec 1 


universal gas constant 
molecular weight of air 
atmospheric pressure 
atmospheric pressure at sea level 
radius of earth 
atmospheric temperature 
log(P/P Q ) 

geometric altitude 
atmospheric density 
earth latitude 


8.311+xiO 3 J/°K 


28.9 6k kg 
N/m 2 
N/m 2 
m 
°K 


m 

kg/m 3 

degrees 


Governing equations 

Atmospheric temperature, density and pressure are related by the 
equation of state 


MP 

p m S 


(A.l) 


and the hydrostatic equation 


dP = ~gp dz 


(A. 2) 


The molecular weight of air M can be considered constant, for the 
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present purposes , at altitudes "below 90 km. 

The acceleration due to gravity g , a function of altitude 
is given "by 



(R o + «)■ 


z 


> 


(A. 3) 


where g , the value of g at sea level, is 

o 

s 9.780356 £l + .0052885 sin 2 4> -.0000059 sin^ (2<J>) J (A. 4) 

and R q , the earth’s radius, is 

R = 6356798 (.9933070 + .0066930 sin <J>P (A. 5) 

o 

Equations (A. 3) - (A. 5) are actually approximations, but they give, 
sufficient accuracy for the model calculations in this study. Equations 
(A. 3) and (A. 4) are those used for the calculations in the 1966 Standard 
Atmosphere Supplements, [2], The earth's radius R Q (<P) at latitude 
, (A. 5), is based on the assumption that the earth 1b an ellipsoid 
with an equatorial radius R q (0) = 6378178 m, and a polar radius R Q (90) 
= 6356798 m. 

Solving for P and p given a temperature-altitude profile 

Let the given temperature-altitude profile be of the form 

|(T lt z i ): z i = 0, 1 90 km | (A. 6 ) 

and assume that for z between z^ and z 


, T(z) is found by 
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linear interpolation. That is. 


T(z) ■ + (z - z ± ) AT^/Az 


(A.7) 


where AT. « T_ , - T, and Az * z. - z. . (Az does not have a 
i i+1 i 1 +1 l 

subscript since Az « 1 km for all i , i - 0, . 90). Assume further 

that between z. and z. _ g is constant and is given by (A. 3) 

i i+l 

with z = z, . 

i 

In order to solve (A.l) and (A. 2) for pressures {P^ and densities 
{p i > it is necessary to have a value of pressure or density 

at some boundary point z^ . Assume that z^ ~ z^ for some i between 
0 and 90 , and that is known at that altitude. 

Substituting for p from (A.l), equation (A, 2) becomes 


dP _ gM 
P KT 


dz 


(A.8) 


If P i is known at then, integrating from z ^ to z i+1 , one gets 


p i+ i - p i 


g(z,)M r i+l 


V? f 

K l 


dz 

T^z) 


(A. 9) 


Using T(z) from (A.7), the integral in equation (A. 9) can be integrated 
to give 

r 

T, 


P s p 
r i+l i 


i+l 


g(z A )MAz 


KAT, 


if A^ f 0 , or 


(A. 10) 
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i+1 


[ g(z i ) MAz 
" KT^ 


(A. 11) 


if AT i = 0 . 

Once the pressures {P.^} are obtained, the equation of state is used 
to get densities. 

MP. 


'i KT, 


(A. 12) 


Solving for z and 0 given a temperature-pressure profile 
Let the temperature-pressure profile he of the form 

| (Tj, Uj) : i = 0, 1, .... 90 | ( A - 13) 

where U = log(P /P ) and P is sea-level pressure. P Q can 

«] V v 

actually he the pressure at any boundary point z^ , not necessarily 

z = 0 . Values of U 4 ranging from U ■ 0 to « -6.U cover 

h J u y 

the atmosphere up to approximately 90 tan. 

It is assumed that for any U betv« 
can he obtained by linear interpolation. That is. 


It is assu m ed that for any U between Uj and 


T(U) * + ( u - Uj) ATj/AUj 


(A.lU) 


where AT^ = T J+1 - ’I J and AU^ = U J+1 - Uj . 
Equation (A. 8) can now be written 


M 


dz 


(A. 15) 


(R o ♦ .)' 
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Integrating from U q 88 0 to , (z^ ■ 0 to z^) * one gets 


K 

M 


i: 


T(u) du = - 


g R z 
o o j 

Tr + z7) 

O J 


(A. 16) 


and solving for z yields 0 
J K 


jl r 

g Ml T(u 

0 ./u. 


) du 


K 


6 MR 
o o 


V 

1 


T(u) du 


(A. IT ) 


where the integral is approximated by 


0 

1 


J-l (T. +1 + T ) 

T(u) du = E | Ail. | o 

i=0 


(A. 18) 


The 


density p at z is again found by using the equation of state. 

<3 0 



APPENDIX B. EQUATIONS FOR CALCULATING SOLAR HEATING RATES 


List of Symbols 


Symbol 

A 

a 

B 

C P 

g 

S o 

K 

k. 

i 

14 

M0 3 

P 

P 

o 



Definition 

shaping parameter in 0^ model 

surface albedo = reflected radiation/ 
incident radiation 

parameter in 0^ model 

specific heat at constant pressure 

acceleration due to gravity 

acceleration due to gravity at sea 
level 

universal gas constant 

linear absorption coefficient for 
step i 

molecular weight of air 

molecular weight of ozone 

atmospheric pressure 

atmospheric pressure at sea level 

partial pressure of ozone 

maximum P0^ (model parameter) 

earth* s radius 

flux of solar energy 

flux of solar energy in step i 

portion of solar constant in step i 

contribution of Rayleigh scattering 
to planetary albedo 


Units /Const ants 


100U.T J/kg/°K 

m/sec 2 

m/sec 

8.31**3xl0 3 J/°K 

-1 

m 

28.96U kg 

1+7.9982 kg 

N/m 2 

N/m 2 

N/m 2 

N/m 2 

m 

W/m 2 

W/m 2 

W/m 2 
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Symbol 

T 

T 

s 

U 

u* 


z 

e 

K 

X 

P 


Definition Units /Constants 

°tr 

atmospheric temperature 

solar heating rate of ozone °K/sec 

height parameter = log(P/P o ) 

# 

value of U where PO^ a PO^ 

(model parameter) 

altitude m 


cosine of the sun* 8 zenith angle 

mass absorption coefficient for step i 

wavelength 

atmospheric density 

ozone density 


m^/kg 

micron p * 10 ^ m 

kg/m 3 

kg/m 3 


Ozone absorption coefficients 

are three spectral bands in which ozone absorbs solar radiation* 

These are the Hartley band with .2 < X < .3 V , the Huggins band with 

*3 < X < .37 \i , and the Chapius band with < X < *7 U • 

The present study uses an absorption coefficient model, [36], 

in which these spectral bands are divided into three "steps" each having 

a constant mass absorption coefficinet ,i = l, 2, 3* The three 

steps and their respective mass absorption coefficients are listed in 

Table B.l . The K.*s given here are based on ratios of step heights 

K : |< 2 : = 1: -02: *0005 calculated from data in Inn and Tanaka, [1*5]- 

The value of was selected to give:. a linear absorption coefficient 

k.(z) , defined by 
i 
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Table B.l 

THREE STEPS IN OZONE ABSORPTION MODEL 




St ep 

Wavelengths , u 

Mass Absorption, 
Coefficient , 

m 2 /kg 

Fraction of 
Solar Constant 11 

HI 




n 

.200 - .300 

1292.0 

.01203 

2 

.300 - .335 

25.84 

.02113 

3 

.335 - .370; .400 - .700 

0.6460 

. 40982 


#Solar Constant = flux of solar energy incident on the "top” of 

the atmosphere =1353 W/m^ 
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k±(z) * (z) (B.l) 

equal to .0005 m -1 at the altitude where ozone density is maximum. 

This value was chosen by Su, [36], to yield local heating rates in good 
agreement with those of an earlier model by Manabe and Strickler, [191 • 
The value of maximum ozone density used for calculating was 

3.87 x 10"^ kg/m 3 as given by the Mid- Latitude Ozone Model in reference 

[ 2 ]. 


Solar heating rates 


The solar heating rate of ozone at altitude z is 




1 dS 
pc p dz 


(B.2) 


where the solar energy flux S(z) is the sum of three components S^(z) 
corresponding to the three steps in the model. These components are 




k^(x) dx 


where the values of S . are based on Table B.l; s_. and a are 

00 , 1 ** 

respectively, 0.07 and 0.10, as used by Su^ [36], and 0 — .5 is the 
cosine of the effective mean zenith angle. 

Taking the derivative of (B.3) with respect to z , one gets 
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dS. \ 

- s »,i (l - S R ) k i (z) exp 


* — ■ j"- g- j" ^(x) dxj 


+ a0 


* [• y £ V*’ **] " p [- ^{l *‘ w “]J , 


Letting k^z) = k^U) and summing (B.U) over all values of 1 , 
gets 


one 


(1- sj p-(z) 3 


R 3 

'F i*l 


= e, p(»r — 2 s - 


,i K i { exp [- T f z p 3 (x) dX ] 


+ a8 exp T- -j- j P 3 (x) dxj exp K i ^ J Q p 3 ^ dx j | 

° (B.5) 

The present ncdel uses U = log(P/P Q ) as its independent variable, 
and, therefore, (B.5) must be expressed as a function of U . The 
equation of state for a constituent density is 


MO^ 

KI 


(B.6) 


and, therefore, the mixing ratio of 0^ is 


P 3 _M0 1 P0 1 
0 MP 


Then, utilizing (B.7) and the hydrostatic equation, one can vrite 


(B.T) 


MO- 1*2 


I P3W dx • - I P0 3 (u) du 

Jz 1 ^ Jv 1 


(B.8) 
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The acceleration due to gravity g , which actually varies slightly 
with altitude (i.e., by about 3 percent between sea level and 90 km), is 
replaced by a constant 

* = g o (R q + 9o7 


which represents its effective mean between sea level and 90 km. 

Using (B.7), (B. 8) , and (B*9), the heating rate can now be written 

as a function of U . That is. 


(1 - S R ) mo 3 po 3 (u) 


V u) “ c, M r exp (u) 


a0 a 


exp 


j, S %i K t Jexp [- tJC/ 0 3 (") dU ] 

J*- jJ P0 3 (u) duj exp Sj P0 3 (u) duj 


where 


K MO 

* i~ 

1 gM 


(B.10) 

(B.ll) 


The use of a constant g “ g , and a constant sea-level pressure 

♦ 

P in (B.10) introduces a small error in the T (U) profile. However, 
o B 

in the temperature model, T (U) is standardized so as to become a 

random variable with mean equal to zero and standard deviation equal to 

one* Thus, any error which affects the mean and/or standard deviation 

of T (U) does not affect the temperature model, 
s 

Calculating optical thicknesses 


The quantity 
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1 

0 



k.(x) 


dx = 



P0 3 (u) du 


(B.12) 


is called the optical thickness or opacity of the ozone "between and 

Z 2 * 

The following calculations are based on the ozone model explained 

in Chapter IV. A value of U located at approximately 50 km altitude, 

denoted , is used as the dividing point for the two different PO^ 

formulas. For U < U c . (i.e., at altitudes above the dividing point), 

50 

PO^ is given by 

P0 3 (U) = exp(AU - AU* + B) (B.13) 


where for a particular profile A , B , and U* are constants. Thus, 


when U < U e< . , the optical thickness above z(U) is 
50 


oo 

1 / 

Jz( U) 


k^x) dx = jg- P0 3 (U) 


(B.l4) 


If U > U_ rt , Green* s formula (Chapter IV) gives 
50 

• e L(u) 

P0 3 (u) = “ P0 3 (1 : jM f 


where 

L(U) = A (U - U*) 


(B.15) 


(B.16) 


50 * 


Thus, for U > U 



CDjM 
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k.(x) dx = ^ 


exp [L(U 50 > + B] 


, UP0 * \ 
A6 


1 + e 




1 + e 


lTuT 


(B.17) 



APPENDIX C, ESTIMATING THE PARAMETERS IN GREEN’S FORMULA 

USING OZONAGRAMS 

List of Symbols 


Symbol 

Definition 

Units /Const ant s 

A 

shaping parameter in Green’s formula 


dq/du 2 

slope of the ozonagram at a point above U* 
where P0^ = q^ 


M 

molecular weight of air 

28.96** kg 

mo 3 

molecular weight of ozone 

1*7.9982 kg 

p 

atmospheric pressure 

N/m 2 

p 

o 

atmospheric pressure at sea level 

N/m 2 

P* 

# 

atmospheric pressure where P0^ ■ P0^ 

N/m 2 

P0 3 

partial pressure of ozone 

N/m 2 

« 

po 3 

maximum partial pressure of ozone 

N/m 2 

% 

ozonagram value of PO^ at height u^ 

N/m 2 


arbitrarily selected value of P0_ where 
slope dq./du 2 is estimated 

N/m 2 

to 3 

total ozone 

atm-cm 

u 

height parameter - log(P/P Q ) 


”l 

arbitrarily selected value of U above U* 
where P0^ = is measured 


u* 

value of U where P = P* 


P 3 

ozone density 

kg/m 3 
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Criteria for a good fit 


In fitting Green 1 s formula. 




(1 + e x ) 2 


(C.l) 


where x = A log(P/P#) , to an ozonagram, ideally one could read the 
coordinates (PC> 3 » P*) off the curve at the location where PO^ is 
maximum. Unfortunately, however, the data usually have an appreciable 

amount of scatter at that point, and, consequently, it is sometimes 

« 

difficult to determine values of PO^ and P* which produce a good 
fit. 

In the present study, a "good” fit of Green’s formula is one which 

(1) fits the smooth segment of the PO^ profile above its maximum and 

(2) has the correct total ozone TO^ as reported on each ozonagram. 


Calculating the shaping parameter A 
« 

If PO^ is known, the shaping parameter A can be selected to 
satisfy this latter criterion regarding total ozone. Total ozone, 
in atmosphere-centimeters, is given by 


to 3 


22U0 

M0 3 


00 



P 3 (x) dx 


(C.2) 


The ozone density p 3 can be expressed in terms of P0 3 to give 


to 3 * 


C 

_ 22^0 C 

5m L 


P0 3 (u) du 


(C . 3 ) 
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where U * log(P/P Q ) . (See equations (B.6) - (B.9) in Appendix B). 
Using Green’s formula for PO^ > one gets 


8960 PO 


TO, 


1 [l - (1 + e^V 1 ] 


(c.U) 


gMA 


and since P » P* (i.e., U* is a large negative number ) , then 
o 

8960 PO, 

TO, = — — * (C.5) 

3 gMA 


This approximate expression is used to calculate A for any PO^ an< ^- 


T0 3 . That is, 


A = 


8960 PO, 

gM TO, 


(C.6) 


(The A q mentioned in Chapter IV is 8960/gM) . 


Calculating PO, and P 


Since the portion of the measured PO^ curve above the maximum 
is generally smooth, points on this portion can be read with little 
difficulty. The following information, read from each ozonagram, is 

ft # 

used to calculate PO^ and P : 

ft . * 

(a) a point (q^,^) on the curve above (PO^jP ) , where 
PO 3 = q 1 and log(P/P o ) = u^ , and 

(b) the slope dq/du 2 - APO^/hU ..at a second point PO^ = q^ 

£ 

above (PO ,’P. ) . 

From thi s information , i . e . , ( q x , 1 ^ ) and ( q g ,dq/du 2 ) , values 

* ft 

^ Q.TQ calculated such that the fitted curve passes 
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through and lias a slope equal to dq/du 2 at the point above 

the maximum where PO^ = . 

The derivative of Green’s formula with respect to U is 


dPO^ A PO 2 (1 - e 2x ) 

du , x 
4P0^ e 


(C.T) 


Substituting for A from equation (C.6) and letting 
and P0 3 = q 2 , (C.T) becomes 



22U0 ^ (1 - e 2x 2) 
gM T0 3 e X 2 


dFO^/dU * dq/du 2 


(C.8) 


Solving for 



one gets 



* 


= - y +%y + i 


where 


gM TO- dq/du 2 
v » . * 1 ■ ■ — 

UU80 q| 


(C.9) 


(C.10) 


A p 

Thus, the value of e X at PO^ * q 2 , e , is a function of known 
quantities and can be calculated. Then, substituting this into Green's 
formula for , it is possible to solve for 



qg(l + e ^) 2 

4 e x 2 


(C.ll) 


where e 2 is given be (C.9)« 

Now, a value for A can be calculated from (C.6), and P* = 


U* 
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can be found by solving Green's formula at (q^,u^) , 


* e 


q, * UPO Q —z 

1 3 (1 + e X l) 2 


for U* = “ x i/a ' (C.12) for e gives 


^ = [vf - /r - l]‘ 


where r = 


PO^/q^ . Thus, 


P* 


p e 1 {Jr - /r - l) 


-2/A 


(C.12) 


(C.13) 


(C.lM 



APPENDIX D. EQUATIONS FOR FITTING PEARSON DISTRIBUTIONS TO DATA 

Let X be a continuous random variable and let f(x) be its 
probability density function. Karl Pearson, [^2], claimed that the 
probability density functions of most statistical distributions arising 
in practice belong to the family of density functions satisfying the 
differential equation 


|f = (x - a) f. (ij.1) 

x b 0 + V + V 

where a , b Q , , and b 2 are population parameters (constants). 

Pearson identified thirteen distribution types which satisfy 
(D.l), These consist ..of the normal distribution pluB twelve types 
which are designated Type I - Type XII. Of these, the normal distri- 
bution, and the Type I, Type III, Type IV, and Type VI distributions 
are considered the major Pearson distributions, whereas the other 
eight are called transitioned distributions, [U6], and can usually 
be treated as special degenerate cases of the five major Pearson dis- 
tributions. Methods are presented here for fitting the five major 
distributions to empirical distributions, (data) . In instances where 
a transitional distribution is identified, it is assigned to one of the 
five major distribution types. 

/ 

Any probability density function of the Pearson family can be 
completely determined by its f6ur moments ,i=l, . .., U, 
defined by 
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00 

■/ 


f(x) dx 


(D.2) 


Assume that a random sample of size n is drawn from the X 
population, and let the population moments he estimated hy the 


sample moments as defined hy 


1 i 

M = - £ X, 

i n J 


(D.3) 


Let M| denote the i central moment in the sample as defined hy 


n 


M! = i Z (X, - X) 1 
1 n j=l i 


(D A) 


2 

where X ** is the sample mean. Then ■ 0 , and - s 
is the sample variance. 

The major Pearson distributions are characterized hy three 
parameters 0^. , $ 2 , and k which are defined as follows: 


3. * 


< M 3 )2 

(M 2 )3 


M i 


(M -) 2 


and 


k = 


(B 2 + 3)' 


1*(46 2 - - 3B X - 6) 


(D.5) 


(D.6) 


(D.T) 


2 

If the sample variance s # 0 , then, 3^ > 0 , $ 2 > 0 , 
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k = Q if and. only if ■ 0 . 

In most cases, the values of 3^ » Bg » K determine the 

Pearson distribution to which the sample belongs. A fourth parameter 
is sometimes required to distinguish the transitional distributions. Figure 
D.l illustrates the various regions in the 3-^ - 3^ pi® 116 which are 

associated with the first seven Pearson distributions. The normal dis- 
tribution corresponds to the point labeled "N" where 3^ * 0 and Bg ~ 3 * 
These first seven Pearson distributions plus the normal distribution are 
listed below with the criteria for their selection and equations for 
estimating their parameters. 

There are two circumstances in which 3-^ and are inadmissible- 

and cannot be used to estimate Pearson distributions. The first is the 
violation of a constraint on all frequency functions which requires 
that 

g 2 - B 1 > 0 (D.8) 

The second circumstance is the violation of a constraint necessary 
in order for 3 1 and 3 g to lead to a Pearson distribution. This 
constraint requires that 

15 B-l - 8 B 2 + 36 > 0 (D.9) 

Three instances arose in the Monte Carlo Shuttle entry study in 
which constraint (D.9) was violated. In these cases an attempt was 
made to find the "closest" Pearson distribution fitting the data. 

To accomplish this, a line L was drawn through the inadmissible point 
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(B , 3 ) perpendicular to the line 15 3-, - 3 3g + 36 = 0 . ^ ee 

Figure D.l) . A point on the line L Just inside the admissible region 

# * 

vas then selected, and its coordinates (3^ > 3g) were used as new 
adjusted values for & and B 2 . The third and fourth central moments , 
and Mj^ , were then adjusted to be consistent with the definitions 
of & 1 and $ 2 in equations (D.5) and (D.6), respectively. 

In the case of one parameter (i.e. , maximum dynamic pressure during 
the summer entries) this method succeeded in producing a reasonable 
Pearson fit for the empirical distribution. In the two other cases of 
inadmissible 3 i 's, (i.e., summer down range and cross range), this 
adjustment produced J-shaped Type VI distributions with poor resemblances 
to the parameter histograms. 

The following summaries explain the criteria used for selecting 
each Pearson distribution and the methods used in estimating their 
parameters. Ho effort is made to show derivations of the formulas. 

If the reader is interested, a good explanation can be found in Kendall 
and Stuart, Vol. 1, [U2]. 


The Pearson Type I Distribution 

(A) Criterion for selection: K < 0 

(B) Probability density function: 

£ > - 1 ; * < 1 
a 2 

elsewhere 


|y (1 + T y 1 (1 - t ) 


X 


f(x) =' 


“1 


(D.10) 
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where x = X - . 

(C) Equations for estimating parameters: 
Let 


and 


The estimates of 


R = 6(e 2 - - 1 )/(6 + 38 1 - 28 2 ) 

(D.ll) 

S = [8 X (R + 2) 2 + l6(R + l)] 1/2 

(D.12 ) 

T = R(R + 2) 8^ 2 /S 

(D.13) 

, m 2 , a x , a g , and Y q are, 

respectively. 

( (R - 2 + T)/2 if Ki < 0 


fiL = < 3 

(D.lU) 

X l (R - 2 - T)/2 if > 0 


*2 = R “ 2 - \ 

(D.15) 

=.5 s S(m 1 + l)/(m^ + ^ + 2) 

(D.16) 

a 2 =.5 b S(fi 2 + l)/( + £ 2 ♦ 2) 

(D.17) 


and 


A 


Y 

o 


r(m 1 + + 2) 

r(m 1 + l)r{m 2 + 1) 







(D,18) 


The Pearson Type II Distribution 


(A) Criterion for selection: k « 0, P 2 < 3 

This is one of the transitional distributions, and is treated as 
a subclass of the Type I distribution. 



ikx 


The Pearson Type XII Distribution 

(A) Criterion for selection: K = 

(B) Probability density function: 


f(x) - 


Y o (1 + l )a-1 ex P(- bx > 


— > -1 
a 


elsewhere 


(C) Equations for estimating parameters: 

The equations for estimating a , b , a , and are, 

a = 

b = 2 s 2 /M^ 


and 


a « a/B 

; _ M a (a • x) exp(-s) 
Y 0 r(&) 


The Pearson Type IV Distribution 

(A) Criterion for selection: 0 < K < 1 

(B) Probability density function: 

f(x) = Y - |X ♦ (f - b 2 exp [ -v t^- 1 (£ - $ ] 

u u r a r 


where x = X - . 

(C) Equations for estimating parameters: 


(D.19) 

respectively , 
(D.20) 
(De2l) 
(D.22) 
(D.23) 


(D.24) 


Let 




R = 6(e 2 - - l)/(28 2 - 3B 1 - 6) 

(D.25) 


S = [l6(R - 1) - - 2) 2 ] 1/2 

(D.26) 

and 

T = R(R - 2) 8^ 2 /S 

(D.27) 


Then the equations for estimating m,» a , v , r , and Y^ are , 
respectively, 


& » (R + 2)/2 


(D.28) 


a “ sS/U 


A 


V 



if M* < 0 
if > 0 


R 


and 


where 


A 


Y 

o 


0.?Q89l*32 f 1/2 ( cos 2 » 1 

rr TTf + 1 ) a? - Hr 

a(cos 9 ) l 



4> 


tan 


-1 



(D.29) 

(D.30) 

(D.31) 

(D.32) 
(D. 33) 


The Pearson Type V Distribution 

(A) Criterion for selection: k = 1 

This is one of the transitional distributions and is treated as a 
subclass of the Type VI distribution. 
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The Pearson Type VI Distribution 

(A) Criterion for selection: K > 1 

(B) Probability density function: 

Y (1 + - y 1 (X + £ ) m£ f > -1 ! f > -1 

w , ° *1 a 2 a l 2 

f(x) = 

0 elsewhere 


where x = X - . 

(C) Equations for estimating parameters: 
Let 


R = 6(6 2 - - l)/(6 + 3^ - 2B 2 ) 

(D.35) 

S = [0 1 (R + 2) 2 + 16 (R + 1) ] 1/2 

(D.36) 

1/2 

and T = R(R + 2) 0^ /S 

(D.37) 

Then the equations for estimating m^ » m^ « a l ' a 2 * 
respectively: 

and Y are , 
0 

^ = (T - R + 2)/2 

(D.38) 

m 0 - (T + R - 2 )/2 

(D.39) 


- D/C^ - fi 2 - 2) 

if 

M 3< 0 

- l)/(m 1 - fi 2 - 2) 

if 

> 0 

i S/2 

if 

< 0 


a. - s S/2 


if > 0 


(D .1*1) 



lUU 


and. 


Y = 


r(m 1 ) 


r(m 2 + - fi 2 - i) 


(f )'W |ti ,A |Si 


m. 


(D.42) 


The Pearson Type VII Distribution 

(A) Criterion for selection: tc * 0 , 3^ > 3 

This is one of the transitional distributions and is treated as a 
subclass of the Type IV distribution. 


The Normal Distribution 

(A) Criterion for selection: 3-^ = 0 3^ = 3 

(B) Probability density function: 

f(x) = — ^ — exp (-x^/2a 2 ) (D.^3) 

0 ^ 2 u 


where x = X - . 

(C) Equations for estimating parameters . 

An estimate of a is s , the sample standard deviation. 
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